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Abstract. Per odic Schur process is a generalization of the Schur process intro- 
duced in [ORl] (math. CO/010705^ . We compute its correlation functions and their 
bulk scaling limits, and discuss several applications including asymptotic analysis of 
uniform measures on cylindric partitions, time-dependent extensions of the discrete 
sine k ernel, and bulk l imit behavior of certain measures on partitions introduced in 
[NO] ( ^icp-th/030623^ in connection with supersymmetric gauge theories. 
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Introduction 

One way to see how the content of this paper is different from previous works 
on the subject is to examine the following three pictures. 
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The leftmost one is a schematic image of an ordinary partition — a way of rep- 
resenting a natural number as an unordered sum of natural summands.^ Partitions 
can also be viewed as ways of inscribing weakly decreasing nonnegative integers 
into unit intervals filling the half-line so that the total number of nonzero entries is 
finite. 

The middle picture represents a plane partition — a way of filling the boxes of 
the square grid in the quarter plane with nonnegative integers so that the numbers 
do not increase as we move to infinity in the directions of the axes, and the total 
number of nonzero entries is also finite. 

Finally, the rightmost picture represents a cylindric partition — a way of filling 
the boxes of the square grid wrapped around a half-cylinder with nonnegative inte- 
gers so that the numbers do not increase as we move away from the border of the 
cylinder in either of the two perpendicular directions of the grid lines. The total 
number of nonzero entries is also required to be finite. 

In this paper we initiate the study of random cylindric partitions. 

Random (ordinary) partitions or, in other words, various probability measures 
on partitions, have been extensively studied since 1940's. The number of references 
is so large that we will not even attempt to list them. An excellent survey of uses 
of random partitions is available in [02]. 

Random plane partitions are less common, partly because they are much harder 
to study. Substantial progress in understanding the uniform measure on plane 
partitions with given norm (=sum of filling numbers) was achieved only recently, see 
[CK], [ORl], [0R2]. In particular, the authors of [ORl] introduced new techniques 
which allowed them to derive determinantal formulas for the correlation functions 
of random plane partitions with weights proportional to q^°™^ < q < \? 

The main object of [ORl] called the Schur process is a generalization of an earlier 
concept called the Schur measure introduced in [01] to handle certain measures on 
(ordinary) partitions. The range of applications of Schur measures and Schur pro- 
cesses, apart from uniform measures on plane partitions, is remarkably broad; exam- 
ples include harmonic analysis on the infinite symmetric group [01, §2.1.4], [BOl], 
Szego-type formulas for Toeplitz determinants [BOk], relative Gromov-Witten the- 
ory of C* [OP], random domino tilings of the Aztec diamond [J2], discrete and 
continuous polynuclear growth processes in one space and one time dimensions 
[PS], [Jl], etc. 

In this paper we introduce and study a generalization of the Schur process which 
we call the periodic Schur process. We derive explicit formulas for the correlation 
functions of this new process and use them to compute the correlation functions 
of the random cylindric partitions with weights proportional to g"o''™. This result 
allows us to obtain various (bulk) limits of these correlation functions as q ^ 1. 
The limiting cases differ by how fast the radius of the cylinder grows comparing to 
I lng|~^, and the results depend in a nontrivial way on the angle between the grid 
lines and the axis of the cylinder. 

We also present two other applications of the periodic Schur process. 

First, we use it to construct an infinite-dimensional family of determinantal point 
processes on which extend the well-known one-dimensional discrete sine process. 

^In this example the number 19 is represented as 19=5-f-4+4+3-|-l+l+l. 

^For such measures the conditional distribution of the plane partitions with fixed norm is 
always uniform (and independent of q). For this reason these measures are often also called 
uniform. 
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This family includes two previously obtained in [ORl] and [B02] extensions as 
special cases. Such abundance of two-dimensional extensions is rather unexpected: 
In all previously known examples, probabilistic models yielded only one extension 
per model, and for most one- dimensional determinantal point processes no more 
than one extension is known. 

The second application is a computation of the correlation functions and their 
bulk scaling limits for a measure on (ordinary) partitions introduced in [NO] in 
connection with certain supersymmetric gauge theories. 

Let us describe our results in more detail. 

The periodic Schur process depends on a natural number N (the period), a 
parameter t, \t\ < 1, and 2N specializations a[l], . . . , a[N],b[N] of the algebra 
A of symmetric functions.^ The process lives on periodic sequences of 2A'' partitions 

A(^) = A(°) D M^'^ C A(i) D • • • C A(^-i) D C A(^) = A(°), 

and it assigns to such a sequence the weight proportional to 

Here Sa/^'s are the skew Schur functions. The proportionality coefficient (which 
we explicitly compute) is chosen so that the sum of all weights is equal to 1. The 
weights can be viewed cither as complex numbers or as formal series in A®^-'^[<;]. 

When t = the partition A*^"^ = A*^^) must be empty in order for the sequence to 
have a nonzero weight, and the periodic Schur process turns into the conventional 
Schur process of [ORl]. On the other hand, if all the specializations a[/c],6[fc] are 
trivial, the periodic Schur process turns into the uniform measure on partitions, 
which assigns to a partition A the weight proportional to f'"^'. 

Denote by Z' the set Z+i. To any sequence of partitions as above it is convenient 
to associate a point configuration (subset) in {!,..., N} x Z' = Z' U . . . U Z' given 

This set determines the sequence (A^^^ . . . , A^^^) uniquely. 

Correlation functions are defined as probabilities, with respect to the periodic 
Schur process,"' that this random point configuration contains a fixed finite set of 
points: 

Pn{Tl,Xi,...,Tn,Xn) = Pvob ^Xj € {X'f'^ " * + ^}i>l I J = 1> • • • > • 

It turns out that the algebraic structure of the correlation functions substantially 

simplifies if one considers a modification of the periodic Schur process which wc call 
the shift-mixed periodic Schur process. It can be viewed as the product measure of 
the periodic Schur process and a measure on Z given by 

ProblS} = const -z^t^ , S eZ, z eC, 

specialization of A is an algebra homomorphism of A to C. 
*One should remember that these "probabilities" do not have to be positive and, moreover, 
may be just formal series of symmetric functions. 
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mapped to the space of point configurations in {1, . . . , N} x Z' via 

In other words, aU points of the random point configuration of the periodic Schur 
process are shifted by an independent integral valued random variable S distributed 
as above. Here ^ is a new complex parameter. 

The normalization constant of the iS-distribution is the inverse of one of the 
Jacobi theta-functions 

^3(^;i) = E^^** = 11(1-*") n {l + t^z){l + tyz), 

SeZ n>l n=i,|,|,... 

see e.g. [Er, 13.19(16)]. We assume that z ^ . . . , so that 6*3(0;*) ^ . 

The correlation functions of the shift-mixed process are defined in the same way 
as those of the initial process, and we denote these new functions as p^^^*. It is not 
hard to see that pn is equal to the constant term in z of ^3(2;; t) p^''*. 

In order to state our first result, we need to introduce more notation. Set 
(im[k] ■= ^Pm{a[k]), bm[k] := ^Pm{b[k]), where p^s are the Newton power sums. 
For any t = 1, . . . ,N denote 

F{r,0=e.pY.ir^i:bn[k] + ^ E bulk] 

rt>l \ fe=l fe=T+l 

k=l k=T+l / 

Theorem A. The correlation functions of the shift-mixed periodic Schur process 
have determinantal form: For any n> 1 and (ri, xi), . . . , {tu, Xn) G {Ij • • • , x Z' 
we have 

P?'^*(n, a;i; . . . ; T„, x„) = det[K{Ti,Xi\ Tj,Xj)]lj^-^ 
where the generating series of the correlation kernel K{a, x; r, y) has the form 



This statement can be understood in two different ways — as a formal identity 
of series in A®^^[t] or as a nmueric equality under suitable convergence conditions. 

Ramanujan's summation formula for I'i/'i-series shows that the two series in the 
formula above are expansions of one and the same holomorphic function in two 
disjoint annuli. This function can be expressed in terms of Jacobi theta-functions, 
see Remark 2.4 below for details. 

For * = Theorem A yields the determinantal formula for the correlation func- 
tions of the conventional Schur process initially proved in [ORl], see also [Jl] and 
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[BR] for other proofs. Note that the argument presented in this paper provides an 
independent proof of this result. 

Let us point out that Theorem A is not particularly trivial even in the simplest 
case of the uniform measure on partitions, which arises when all specializations a[k] 
and b[k] are trivial. Then one readily sees that off-diagonal values of the correlation 
kernel vanish, and the statement reduces to the fact that the shift-mixed version of 
the uniform measure on partitions is equivalent to the product of countably many 
independent Bernoulli measures, see Corollary 2.6 below. The author was not able 
to locate this fact in the literature, although certain formulas equivalent to it can 
be found in [01]. 

The correlation functions of the initial periodic Schur process are not determi- 
nantal. Nevertheless, they possess a nice multivariate integral representation given 
in Corollary 2.8 below, apart from the fact that they can be extracted from the 
formula of Theorem A by taking the constant term in z as mentioned above. 

The proof of Theorem A that we present in this paper is a verification rather than 
a derivation of the formula for the correlation functions. The initial proof involved 
the formalism of the Fock space and was similar in spirit to the derivations given 
in [01] and [ORl]. However, we decided to leave it out of this paper because of its 
length and c;crtain technical difficulties in justification of formal manipulations with 
operators in Fock spaces. As a matter of fact, our initial inspiration came from the 
work [Ts], where the universal characters - analogs of the Schur symmetric func- 
tions for nonpolynomial representations of the unitary groups were represented 
as matrix elements of certain operators in Fock spaces. We hope to return to the 
Fock space formalism in a subsequent publication. 

Our second result describes the "bulk limit" of the correlation functions of the 
periodic Schur process and its shift-mixed version as t ^ 1 and the period A'', as 
well as the specializations a[k] and b[k], remain fixed. 

Theorem B. Assume that z ^ M<o; aTO[fc]5''m[fc] = 0{BJ") as m — > oo for some 
Q < R<1 and all k = 1, . . . , N; and 



Then as t 1, the correlation functions of the periodic Schur process and its shift- 
mixed version have a limit in the following sense: Choose Xi(t), . . . , x„(t) G 1! such 
that \ \nt\-Xk{t) 7 for allk = 1, . . . , n and some 7 G R, and all pairwise distances 
independent of t. Then for any I < ti, . . . ,Tn < N 



N 



N 




m 



limpff*(Ti,Xi(t) 



;...,T„,a;„(i))=det[4:;?)(: 




=1' 




where the limit correlation kernel has the form 
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Note that the Umit correlation functions are invariant with respect to simulta- 
neous shifts of the "space variables" Xi. 

The assumption Am = -Bm is crucial here; without this condition being satisfied 
even the integrals above may become meaningless because the denominators would 
be allowed to vanish on the integration contours. 

It is rather unusual that the limit can be computed in such generality. 

Before proceeding to cylindric partitions, let us describe in more detail the two 
other applications of Theorems A and B mentioned earlier. 

Recall that a classical theorem proved independently by Aissen Edrei Schoen- 
berg- Whitney in 1951 [AESW], [Ed], and by Thoma in 1964 [Th], states that a 
sequence {c„}^oi cq = 1, is totally positive^ if and only if its generating series has 
the form 

for certain nonnegative parameters {A} and 7 such that ^^^(ai + A) < 00. 

Corollary 1. For any doubly infinite sequences of totally positive parameter sets 
{{oS-''\ li^^\^^^'')}k& O'l^d {{a^''\ l3^''\^^''^)}kez satisfying the additional conditions 

and any c G (0,7r), there exists a determinantal point process^ on Z x Z with the 
correlation kernel 

1 T ^ \-i d( 



IC{a,x;T,y) = 



k—r+1 S> 



where the first formula is used for a < t, the second formula is used for a > t, and 
both integrals are taken over positively oriented arches of the unit circle. 

The equal time values of the kernel above are exactly those of the discrete sine 
kernel on Z: For any r € Z 

1 /"^ dC sin(c(a; — v)) „ 
K.(T,x;T,y) = / — ^ = — r^, x,yGZ. 

Thus, the kernels JC{(7,x:T,y) are extensions of the discrete sine kernel. 

(k) 

The choice of = 1 and all other parameters being zero brings us to the 
incomplete beta kernel of [ORl]. On the other hand, taking 7^*^^ = 7^'^' with all 
other parameters being zero yields the extension of the discrete sine kernel obtained 
in [B02, Theorem 4.2]. 



^By definition, this sequence is totally positive if all minors of the matrix [ci—j]ijy.o are 
nonnegative. 

®This time all the probabilities are nonnegative! 
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As for the second application, wc consider a probability measure on the set of 
all partitions given by the formula, see [NO, §6.2], 



Here /x € and t G (0, 1) arc the parameters, the product is taken over all 
boxes of the Young diagram A, and h{0) denotes the length of the hook rooted at 
the box □. 

One remarkable feature of this measure is that it interpolates between the uni- 
form measure on partitions, which appears at /i^ = 0, and the (poissonized) Plancherel 
measure on partitions (see, e.g., [BOO]), which is obtained from M^^^ by the limit 
transition t 0, iJ ^ ±ioo, ^ ^ > 0. 

The measures Af^ ^ may be viewed as special cases of the periodic Schur process 
with period N = 1, and Theorem B leads to the following statement. 

Corollary 2. As i — > 1. the correlation functions of M^^^t have the following limit: 
Choose Xi{t),.. . ,Xn{t) G Z' as in Theorem B. Then the correlation functions con- 
verge to determinants of the limit correlation kernel 



Now let us return to cylindric partitions. 

It is convenient for us to represent cylindric partitions as periodic sequences of 
ordinary partitions by reading the filling numbers along the diagonal rays which 
form 45° angle with grid lines. For example, the visible part of the cylindric par- 
tition represented by the picture in the beginning of this introduction gives the 

sequence of partitions 



. . . (3) D (1) C (5, 1) D (2) C (7, 2) D (3, 1) C (4, 1) D (1) C (2, 1) D (2) C (7) . . . 



The condition of filling numbers not increasing along the grid lines is equivalent 
to neighboring partitions in such a sequence having interlacing parts. It is also 
equivalent to saying that the Young diagrams of any pair of neighboring partitions 
are different by either adding or removing a horizontal strip. Such a relation be- 
tween two partitions (k, u) is denoted as k y v or k -< u, depending on which of 
these two partitions is larger. 

The choice of ;^'s and ^'s between neighboring partitions is exactly the choice 
of the boundary profile of our cylindric partition near the cut of the cylinder. We 
will fix such a profile by providing two periodic sequences {A[fc]} and {i3[fc]} of O's 
and I's such that A[k] + B[k] = 1; the ;^'s correspond to A[k] = 1 and B[k] = 0, 
and -<'s correspond to A[k] = 0, B[k] = 1. The ratio of the total number of ^'s 
in a period over the total number of )^'s in a period will be called the slope of the 
profile. The slope depends only on the angle between the grid lines and the axis of 
the cylinder. The case >c =1 corresponds to the diagonal rays being parallel to the 
cylinder axis. 

The following picture represents a cylindric partition with slope x <1, and the 
visible part of the boundary profile corresponds to the sequences 




Ae Y. 




{A[fc]} = (..., 1,0, 1,1, 1,1, 1,0,1,...), 
{S[A;]} = (...,0,1,0,0,0,0,0,1,0,...). 
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We will denote by N the period of the sequences above. For any integer m, let 
m{N) be the smallest positive integer such that m = m{N) mod N. 

Proposition. For any profile {A[fc]}^^, O'^ above we have 

1 T-r 1 



E 



cylindric partitions tt 

of the given profile 



nr 



n 



n>l 



1 _ s(p-g)(JV) + (n-l)]V 



pel.JV: A[p] = l 

rjeUV: B[q] = l 



Here \ ■ \ denotes the norm { =the sum of filling numbers) of cylindric partitions. 

We were unable to find this formula in the literature, although it may well fol- 
low from more refined results of [GK] . In the limit when the radius of the cylinder 
becomes large, if = (1, . . . , 1, 0, . . . , 0), B = (0, . . . , 0, 1, . . . , 1), and x remains 
bounded away from and oo (this means that the boundary of the cylindric par- 
titions locally looks like the boundary of the quarter plane), the formula of the 
proposition reproduces the celebrated formula of MacMahon for the sum of the 
weights s""""™ over all plane partitions. 

The applicability of the periodic Schur process to cylindric partitions follows 
from the following basic property of the skew Schur functions evaluated at a single 
indeterminate x: Si^/^{x) = a;!'^'"'''! if k )~ u and otherwise. One readily checks 
that the periodic Schur process with f = and specializations a[k] and b[k] being 
the evaluations at s'^A[fc] and s~''B[k], respectively, for all A; = 1, . . . , TV, is exactly 
the measure const -s"""^™ on cylindric partitions with a fixed profile described by 
{^[fc]} and Prom now on we will use the term "uniform measure" for this 

distribution. 

The above observation implies that Theorem A provides a determinantal formula 
for correlation functions of the shift-mixed modification of the uniform measure on 
cylindric partitions.^ Our next goal is to explain what happens to Theorem B. 

In order to state our next result, we need to introduce a curve Tj^ in the complex 
plane via (here x is a positive parameter) 

I i~:n 1 — ^ — I 

^e^'^|^G[-7r,7r]L 



sm 



l+xr- 



This is a piecewisc smooth closed curve which has a corner-like singularity at the 
point 1. We orient F^^ counterclockwise. 

For a < T set A{a,T] = A[k], and similarly for B{a,T]. The the slope 

X is equal to B{0, N]/A{0, N]. 

^By correlation functions of the uniform measure on cylindric partitions we mean the corre- 
lation functions of the corresponding periodic Schur process, and similarly for the shift-mixed 
versions. 
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Theorem C. In the limit s — > 1, the correlation functions of the uniform measure 
on cylindric partitions with a given profile {A[k], B[k]}^^i have a limit in the follow- 
ing sense: Choose xi{s), . . . ,a;„(s) € Z + i such that as s 1, N\ lns| • Xk{s) 7 
for all k = 1, . . . ,n and some 7 S R, and all pairwise distances Xi — xj are inde- 
pendent of s. Then for any ri, . . . , r„ e {1, . . . , A''} 

^nm_ /?„(ri, a;i(s); . . . , r„, a;„(s)) = det[}CpJ^. {xi - xj)] ^.^^ 



where the correlation kernel has the form 

(1 _ ()'B(<^.-rl(l _ (^-l)^(<^,i-] 



— f - 

p -L 



27ri J 1 
r i 



(1 - C)--^(^''"l (1 - (-i)-Mr,'^] dC 



a < T, 



a > T. 



The function (1 - (J)^(O'-'^] (1 — ( i)^(o.^] takes nonnegative values on T^, and 
thus the integrals are correctly defined. 

The limit density function IC^J{0) does not depend on t, which means that it 
is invariant with respect to rotations of the cylindric partitions. Note also that 
it depends on the profile only through A{0,N] and B{0,N], or, in other words, 
through the period N and the slope >c. 

Interestingly enough, a formal application of Theorem B to the case of cylindric 
partitions may produce an incorrect answer if >f ^ 1. This happens because for 
X ^ 1 the condition Am. = Bm of Theorem B is violated, and we need to deform 
the integration contours in order to perform the asymptotic analysis. 

In Theorem C we kept the period N finite while sending s to 1. The next level 
of difficulty is to consider the case of periods growing together with | lns|~^. 

If the growth of the period is slow in the sense that A''|lns| still tends to 0, 
then we prove that the limit behavior of the correlation functions can be read off 
Theorem C above by taking the limit N 00 and keeping the slope x fixed. In 
the limit one sees extensions of the discrete sine kernel as in Corollary 1 above with 
the parameters a^'^^ or a^*^^ taking values between and 1 and all other parameters 
being zero. Details can be found in §6. 

The case of the period A'' growing in such a way that the product Afj lns| has a 
finite limit, is substantially more complicated. The reason is simple the limiting 
behavior depends on the details of the profile rather than just on its slope. In this 
paper we only consider the case of the slope being equal to 1, and we prove two 
results. 

First, we show that if the profile sequences {^[fc]} and are periodic with 

a finite period (in addition to being periodic with the growing period N) then the 
limit behavior is just the same as in the case of the slowly growing periods. 

Second, we consider the corner-like profiles with {A[A;]} and consisting 
of one block of O's and one block of I's. We compute the limit of the correlation 
functions near two "corners" , where the partitions arc the largest and the smallest. 
The results in both cases are governed by the incomplete beta kernel with the 
density functions given by certain analytic expressions involving elliptic functions. 
Details can be found in §7. 
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To conclude the introduction, let us mention two circles of questions which we 
do not discuss in this paper, but which are certainly very interesting. 

The existence of limit correlation functions which decay fast enough when the 
distance between the arguments grows, is a strong indication for the existence of a 
limit shape of the corresponding random (ordinary, plane or cylindric) partitions. 
Such a decay is present in all the cases we considered. The limit density function 
allows one to predict what the limit shape would look like (see, for example, Com- 
ment 3 after Theorem 3.1 below), but one needs additional arguments to actually 
prove the concentration phenomenon. For example, for the measures M^^t consid- 
ered in Corollary 2 above (and, in fact, for a substantially larger class of measures 
on partitions), the existence of the limit shape as /i — > ±ioo was proved in [NO] by 
variational techniques. 

Also observe that in this paper we consider only what is usually called "the bulk 
scaling limit" of the correlation functions. It would be very interesting to study the 
"edge scaling limit" as well. In the case of the uniform measure on (skew) plane 
partitions different edge scaling limits were computed in [0R2]. 

Acknowledgements. The author is very grateful to Grigori Olshanski for numer- 
ous discussions which were extremely useful and stimulating. In particular, a large 
part of §8 below consists of results of such discussions. The author would also like 
to thank Eric Rains for pointing out a Cauchy type determinant formula for theta 
functions due to Frobenius; it is used in the proof of Corollary 2.8 below. 

This research was partially supported by the NSF grant DMS-0402047 and the 
CRDF grant RIM1-2622-ST-04. 

1. Periodic Schur process 

Fix a natural number N and consider a measure on periodic sequences of 2N 
partitions 



by specifying the weight of such a sequence to be equal to 

W(A, /z) = ls^(o)/^(i) (a[l]) Sa(i)/m<i) (^[1]) • • • Sx(N-i)/^(N) {a[N]) s^(n) /^(n) {b[N]). 

Here Hs a parameter and a[m], b[m], m = 1, . . . , N , are arbitrary specializations of 
the algebra A of symmetric functions.^ We will use the notation 



where pk's are the Newton power sums. 

The weights W(A,/i) may also be viewed as elements of A®^^[f]; in that case 
the notation ,f{a[k]) or ,f{b[k]) for / G A just indicates to which of the copies of A 
in A®^^ this symmetric function belongs. 

The partition function of such a measure will be denoted as (the notation Y 
below stands for the set of all partitions including the empty one) 



= aW d m^'^ C A(i) d • • • C A(^-i) d /x(^) c A(^) = AW 



ak[m] := lpk{a[m]), bk[m] := lpk{b[m]), 






*A specialization of A is an algebra homomorphism of A to C. 
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We will use the term periodic Schur process for this measure. 

It is clear that without loss of generality we may consider only the distribution 
of A's. Indeed, by making some of the specializations a[/c], trivial^ we can force 
any given fj!-'^ to coincide with a neighboring A*^ ' ^ . 

The periodic Schur process as defined above is not symmetric with respect to the 
circular shifts A^^^ i-^ A(''+' '""'^ This lack of symmetry can be easily repaired. 
Introduce new specializations a[m],6[m] of A by assigning the following values to 
the power sums: 

OfcH := s"'''"ofe[m], hk[m] := s'^^ftfeN, 
where =t. Then we obtain 

w(a,m) = sI^'''i+i^'''i+-+i^'"'i 

X sxm/f,(i) (o[l]) Sa(i)/^(i) (^[1]) • • • Sx(N-i)/i^iN) {a[N]) Sa('v)/^(") (b[N]), 

and this expression already possesses the rotational symmetry. The index m of 
a[m] and b[m] can now be viewed as an element of Z/NZ. We will use both this 
form of the measure and the initial non-symmetric one. 

When all the specializations are trivial the measure concentrates on the sequences 
with coinciding terms: all A's and /u's become equal. The resulting distribution on 
one copy of Y is the so-called uniform measure: W(A) = tl'^L This observation 
shows, in particular, that if one sets t = 1 then the partition function may become 
infinite. 

lit = then A^"^ = A^^^ must be empty, and the periodic Schur process coincides 
with the conventional Schur process introduced in [ORl], see also [BR]. 

In general, the periodic Schur process may be viewed as the conventional Schur 
process on Y started at the uniform distribution instead of the empty partition, 
and conditioned to yield periodic trajectories. 

When the period N is equal to 2, the periodic Schur process can be viewed as 
an analog of the Schur measure of [01] for the universal characters of the unitary 
groups (see [Ko] for details on universal characters): 

^ W(AW,A(2);M«,m(^))=5[,u),a(=)'](«[1],M1]-)V),a(=)'](«[2]",M2]). 

/n(i),Ai(2)eY 

Here the superscript "— " in a specialization c stands for the change of signs of all 
the power sums: Pk{c~) := — Pfe(c). 
Introduce the notation 

N N 

M ■■= ^ ak[m], Bk := ^ bk[m\. 

m=l m=l 



^The trivial specialization of A is characterized by the fact that all power sums p^, fc > 1, 
specialize to zero. 
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Proposition 1.1. The partition function of the periodic Schur process has the form 
Z{N,t,a,b) = \{Y^^^v[Y.4 E ""[^KW + ^Yzf^) I 

n>l \n=l \l<l<k<N / j 

(oo N ^ 

ra=l k,l=l 

Remark 1.2. The equality above can be viewed either as an identity of formal 
series in A®^''^[t] or as a numeric equality under the assumption that |t| < 1 and 
the scries 5]]„>i "-an[fc]fcn[^] are absolutely convergent for any k,l = 1, . . . ,N. 

Remark 1.3. Here is a different way to write the formula for the partition function. 

For two specializations a and /3 of A set 

H{a;p) = J2sx{a)sx{P). 

If a and /3 are evaluations of symmetric functions at variables (q:„) and (/3„) then 
the Cauchy identity, see e.g. [Macd, §1 (4.3)], reads H{a;P) = Yliji^ ~ It 
is also not hard to verify that H{a; (3) = exp p„(a)p„(/3)/n. Hence, 

1 ^ 
ZiN,t,a,b) = n n H{a[k];b[l]) H 11 H{tM%b[l]) 

n>l l<l<k<N n>lk,l=l 

= n (Yzh^u^i'^'~'^''^[k]rb[i])]iH{s(^^^^ 

n>l \ k>l k<l 

where the multiplication of a specialization a by a scalar q is defined via 

Pn{q- a) := q"pn{a), n>l. 

Proof of Proposition 1.1. The argument uses the same idea as [Macd, Ex. 1.5.28]. 
It is more convenient to work in terms of s, a, b. We have 

A,// 

X Sa(")/a«(i) (a[l]) Sa(i)^(i) (Ml]) ■ ■ • ifi[N]) Sa(jv)^(«) (b[N\) 

X Sx(N)/^(i){sa[l])sx(i)/^,w(b[l]) ■ ■ ■ Sx(N-i)/^(N){sa[N]) Sx(N)/^(N)(b[N]) 
= H{sa[l];b[N])H{sa[2];b[l]) . . . H{sa[N];b[N - 1]) ^ 

X /k(^) (b[N]) s^(i) (sa[2]) • • • s^(iv)/^(iv-i) (b[N - 1]) s^(iv) (sa[l]) 
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where to sum over A's we used the well-known formula 

X] ■sa/^(")sa//j(/3) = H(a; j3) ^ s^/K(/3)sp/«(a), 

A K 

see [Macd, Ex. 1.5.26]. Applying the same trick to sum over /it's we get 

Z{N,t,a,b) = H{sa[l]:b[N])H{sa[2];b[l]) . . . H {sa[N];b[N ~ I]) 

X H{s^a[l];b[N - l])H{s^a[2];b[N]) . . . H{s^a[N];b[N - 2]) 

(l)|-Ll^(2)|j_...j_|^(iV)| 



X 

X s 



^(N)/p(i) (s a[2]) s„(i)/p(i) ib[N]) ■ ■ ■ s^(N-i)/p(N) {s a[l]) s^m/pW {b[N - 1]). 
Continuing in the same fashion after N summations we obtain 

Z{N,t,a,b) = ]jH{s''-^^%^l])'[[H{s^+''-^a[k]-,b[l]) ■ Z{N,t,s^a,b). 

k>l k<l 

Iterating the above procedure we arrive at the final formula noting that both in the 
formal and analytic settings 

lim Z(N, t, fa, b) = Z(N, t, trivial, 6) = TT • □ 

n>l 



Remark 1.4. Using the formula 

= H{l3; 7) X s^/p{a, j)s^/p{P, S) 
p 

one readily sees that the projection of the periodic Schur process to 

is again a periodic Schur process with a shorter period M and modified specializa- 
tions 

fei k2 N 

0'n[i-] = ^an[i\, an[2]= X a„[i], a„[M] = ^ a„[i], 

i—l i—ki-\-l z— ^Af— 1 + 1 

fei k2 N 

&n[l] =X]6n[i], bn[2]= £ 5„[z], 6„[M] = 6„[i]. 

This fact can be used to define periodic Schur processes with continuous time. 
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2. Correlation functions 

For any A e Y set 

£(A) = {Ai-i + i}.>^cZ' = Z+i. 

By definition, the nth dynamical correlation function of the periodic Schur process 
is the probability that the random point configurations C{X^'^''^) contain some fixed 
points Xk & for all = 1, . . . , n: 

Some of the time moments Tfc may coincide, but if {Ti,Xi) = {Tj,Xj) for i ^ j then 
the correlation function vanishes. 

In order to compute the correlation functions it is convenient to introduce a 
modification of the periodic Schur process which we call the shift-mixed periodic 
Schur process. 

The shift-mixed process is a measure on point configurations in the disjoint union 
of N copies of Z' obtained as follows. Let us take in the ith copy of Z' the point 
configuration C{\^^^), where (A^^\ . . . , A^^^) form the periodic Schur process, and 
let us shift all these N point configurations simultaneously by a random integer S 
distributed according to 

Prob{5} = z^t^ /ei{z;t). 

Here 2; ^ is a new parameter, and 9s{z;t) = "^sez^ ''~ partition 

function of the weights z°t~ which also happens to be one of the Jacobi theta- 

functions, see e.g. [Er,13.19]. 

To summarize, in the shift-mixed periodic Schur process the weight of the point 
configuration of the form 

{5 + A«-. + i}^>^U...u{5 + Af'-^-HiK>, 

is equal to z''^t^ W(A, /i), all point configurations not of this form have weight 
zero, and the partition function is equal 6*3(2; t)Z{N, t, a, 6). 
Since 

e3{z;t)=Y[il-n- n {l + t^z){l + e/z), 

n>l „=1 3 5 

see e.g. [Er, 13.19(16)], we will always assume that z ^ —t^^, —t^^, ■ ■ ■ , so that 
the partition function is never zero. 

The dynamical correlation functions of the shift-mixed process are defined in the 
same way as those of the initial process; we will denote them by p^^^*- 

Proposition 2.1. The dynamical correlation functions of the periodic Schur pro- 
cess and its shift-mixed modification are related as follows: 

pn{Ti,xi; ...■,Tn, x„) = Constant term in 2; of { dz{z\ t) p^'^'^*(ti, xi; . . . ; r„, x„) } , 

"'■^^^^{Ti,Xi;...;Tn,Xn) = V Z'^t"^p„(Ti , Xi - S] . . .]Tn,Xn - S). 



sez 
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Proof. Follows directly from the definition of the shift-mixed process. □ 

The reason for introducing the shift-mixed process is the fact that this process is 
determinantal, i.e., its correlation functions can be written as certain determinants. 
To make an exact statement we need additional notation. 

For any t = 1, . . . ,N set 



N 



n>l 



k=l 



fc=r+l 



k=l k=T+l / 



(1/0" 

nN 



N+T 



1-S 



,kn 



fc=T + l 



Theorem 2.2. The dynamical correlation functions of the shift-mixed periodic 

Schur process have determinantal form 

pT^\ti,xi; . . . ; t„, a;„) = det[K{Ti,Xi; Tj,Xj)\yj^-^ 
where the generating series of the correlation kernel K{a, x; r, y) has the form 



K{a,x;T,y)Cv^ = < 



( F{a,C) ^ (C^)"' 



a > T. 



Remark 2.3. Similarly to Proposition 1.1, the formula above carries two state- 
ments: One holds in the algebra of formal series in A®^^ [t] with z being an arbitrary 
nonzero complex number. To decompose the right-hand side into a series in t one 
uses the expansions 

1 ^ f Ei>o(-^i'")S > 0, 1 ^ ^ 1_ 

l + zt^ y-Y,^^^{-zt"')-\ m<0, l + {zt^)-^ l + zt^' 

The other statement is a numeric equality which we prove under the following 
convergence conditions: 

an[k],hn[l] = OiR^) for some < i? < 1 and all l<k,l<N- 
\t\<l; ZT^O,-t^^,-t^^,.... 

These conditions guarantee that the generating function above is an analytic func- 
tion in ( and r] varying in an annulus either slightly outside or slightly inside the 

unit circle (depending on whether a < t or a > t), and the values of the kernel are 
obtained as the Laurent coefficients of this function; see also the next remark. 
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Remetrk 2.4. The extraction of coefficients of the generating series can be per- 
formed by computing the corresponding contour integrals. Also, both series in 
the formula for the kernel above sum up to one and the same analytic function in 
disjoint domains: 



n>l 



-y/yOi{yz\t) 



= < 



^ T 



y" 



meZ' 



^ T 



ym 



toGZ 



i< |y| < \t\-\ 
1*1 < |y| < 1- 



These are special cases of Ramanujan's summation formula for iV'i-series, see e.g. 
[GR, §5.2]. Indeed, the second sum equals 



^^^ + zt"'+i 1 + zti ■^r, (-zti;t) 



rE 



-iilji{-zt^;-zt^;t,y) 



1 + Zt2 

^ {t,t,~tizy,~ti/{zyy,t)oo _ {t;t)i,y/y e3{yz;t) 



l + zt2 {-t2z,-t2/z,y,t/y;t)oo O^i-yt 2-t)e3{z;t) 
where we use the conventional notation 



{a; t)m = n (1 - '^*")' (^^5 = 11 (1 - 



n=0 



n=0 



The first sum is obtained by the change {y, z) i— > {y~^, z~^). 

Hence, in the analytic setting the formula for the correlation kernel above can 
be rewritten as follows: 



K{a,x;T,y) 



n„>i(i-t"f 

03{z;t) (27ri)2 



F(a,C) 03{zCr,;t) 



dC^drj 



where both integration contours are simple loops going around the origin in pos- 
itive direction such that R < |C|, jryj < R~^, and for a < t we have 1 < \(ri\ < 
min{i?~^, while for a > t wc have max{i?, |t|} < \(rj\ < 1. 

Remark 2.5. In the limit t ^ both the periodic Schur process and its shift- 
mixed version turn into the conventional Schur process of [ORl]. Accordingly, 
Theorem 2.2 yields a dcterminantal formula for the correlation functions of the 
Schur process derived in [ORl] (see also [Jl] and [BR] for different proofs). Note 
that the summation formulas of Remark 2.4 just turn into geometric series 



yy 



= < 



\y\ > 1, 



E 2/", |2/I<l- 
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Corollary 2.6. The shift-mixed uniform measure on partitions is isomorphic to 
the product of independent Bernoulli random variables on {0,1}^ = {{Xm)me2'} 
with 

Prob{a;„ = 1} = ' ^ 

in the following sense: Pairs (A, 5) G Y x Z are in one-to-one correspondence with 
the sequences (xm) € {0, 1}^ of nonzero weights via 

1, if mG {5 + Ai-i+i}, 
0, otherwise, 

and the weight of such a sequence is equal to 
Prob{A,5} = 



Proof. Follows from the fact that for the trivial specializations a and b the periodic 
Schur process turns into the uniform measure (sec §1), the function F(- , • ) becomes 
identically equal to 1, and the equal time values of the correlation kernel are readily 
seen to be equal to 

1 + zf" 



K{T,x;T,y) = , , ■ □ 



Remark 2.7. Corollary 2.6 is equivalent to [01, (3.14)]. It is also fairly easy to 
prove this statement independently by explicitly computing the weight of a sequence 
in {0, 1}^ . This will essentially be done in the proof of Proposition 2.12 below. 

Before proceeding to the proof of Theorem 2.2 let us draw one more corollary. 

Proposition 2.1 explains that one can obtain the correlation fimctions of the periodic 
Schur process by extracting the constant coefficient in z from the determinantal 
formula of Theorem 2.2. In fact, this extraction can be performed explicitly yielding 
a multivariate integral representation. To state the result it is more convenient to 
work with another Jacobi theta function 9i{x;t) defined as follows, cf. [Er, 13.19]: 

oo 

0i{x;t)= J2 i-^Tt^^^x'^+i ={t,t)^{xi -x-i){tx;t)^{t/x;t)^. 

n=—oo 

Since, as was mentioned above, 9s{x;t) = {t]t)oo{—Vix\'t)oo{—Vi/x]t)oo-) we have 

Corollary 2.8. The dynamical correlation functions of the periodic Schur process 
in the analytic setting can be written as 

Pn('^l) Xi ; . . . 5 Tn-Xji^ 



+1 
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where the integration variables Q and rn range over circles \Q\ = ai, \r]i\ = (3i such 
that 

1 1 1 

ai> — >Q;2>-3->•••>Q;„> — 
Pl P2 Pn 

and all the radii and are close enough to 1. 

Proof. We start with the interpretation of Theorem 2.2 given in Remark 2.4 and 

note that 



e3{z;t)0s{-Cvt--^;t) di{-zVi;t)9i{Cr,;t) ei{z;t) 0i{Cr,;t) ' 

where wc used the notation z = —z\fl. The following Cauchy-type determinantal 
formula is due to Frobenius, see [R, Lemma 4.3], [F]: 



det 



i{zCiVj;t) 



Since 



0i{z;t) e3iz;t)^/Cir]i . ..(nVn 

and the constant term in z of ^3(-zn"=i dVi'^i) is equal to 1, the application of 
Proposition 2.1 concludes the proof. □ 

Note that Corollary 2.8 implies the generating series of the density function of 
the uniform measure on partitions is, up to a constant, the inverse of the first Jacobi 
theta-function: 

5^pi(a;)r = (t;i)cx>/ei(e;i)- 

Proof of Theorem 2.2. We will provide a proof for the algebraic variant of the theo- 
rem when both sides are considered as formal series, see Remark 2.3. The numeric 

equality in the analytic setting is a mere corollary: Under the convergence condi- 
tions stated in Remark 2.3 the formal series in both sides are absolutely convergent, 
and hence the fact that they coincide termwise implies that their sums are equal. 

Our proof is based on the following well known statement: 

Let X be a finite set; let L be a \X\ x |X| matrix with rows and columns marked by 
the points of X with matrix elements from an algebra A, and assume that det(l-l-L) 
is an invertible element of A. Consider an >l.-valued measure on the set 2^ of all 
subsets of X given by 

where Lx is the symmetric submatrix of L corresponding to X: 



Lx = \\L{xi,Xj 



\xi.XjGX- 



Then the correlation functions of this measure are also given by minors of a matrix: 
For any Y C X 

p{Y) = Prob{X \ Y CX} = detKy, 
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where K = L{1 + L)-\ 

Proofs of this statement can be found in [Mace, DVJ, BR]. Measures of the form 
above are often called L-ensembles. 

Interesting examples of L-ensembles usually involve an infinite set X, and thus 
the above linear algebraic statement needs to be adjusted to the specific situation 
at hand. Our case is of the same nature. 

Let us take X = Z' U • • • U Z' (total ot N + 1 copies) and consider the matrix L 
which has the following block form corresponding to this splitting of X: 



- 





. 





Q 


-m 





. 











-m 


. 








. 





. 


. -L[N] 






The Z' X Z' matrices Li, . . . , Ljv are Toeplitz and their matrix elements are given 

by 

{L[k])^y = coefficient of C"^ in exp^(a„[fc]C"" + bn[k]c), 

n>l 

while the matrix Q is diagonal: Q^y = zt^Sxy 

As the algebra A we choose the algebra of formal series in A^'^^^^ [t^^] which have 
at most a finite order pole at t = (in other words, the degrees of t entering any 
element of A must be bounded from below) . Wc will also be using the subalgcbra 
Ahoi of A which consists of series holomorphic at t = (i.e., the series which do not 
contain negative powers of t). The algebra Ahoi has a natural Z+ -filtration induced 
by the degrees of symmetric functions and polynomials in t. Wc will denote its 
filtered components by Ahoi{d), d > 0. That is, Ahoi{d) consists of series whose 
terms are of degree at least d. The algebra A also has a natural topology: Two 
series are close if their difference is in Ahoi {d) for d large enough. With respect to 
this topology the algebra A is complete. 

Recall also that the parameter z entering Q is considered numeric, which means 
that it docs not contribute to the degree coimt. 

The connection between the matrix L above and the periodic Schur process is 
explained by the following statement. 

Lemma 2.9. For any partitions A and v, an integer I > max{£{X),i{i/)}, and any 
fc e {1, . . . , N} we have 

det[(L[fc]).._,,.^._,];_.^^ 

= H{a[k];b[k]) J2 Sx/Mk]K/Mk]) + 0{2l - £{X) - + 2), 

where we use the notation 0{d) for elements of Ahoi{d). 

Proof. To simplify the notation we will omit the index "[fc]" in the formulas below. 
By [Macd, Ex. 1.5.26] we have 



H{a; b)^sx/i^{a)s^/^{b) = ^Sp/x{b)Spi^{a). 
n p 
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Let us split the sum in the right-hand side into two parts — over partitions p of 
length < I and > I. The second part has terms of degree at least 2(1 + 1)— i^Xj—iiy) 
and thus can be ignored. The first part can be rewritten using the Jacobi-Trudi 
formula (see [Macd, §1 (5.4)]) as: 

/9l>P2>--->Pl>0 

where /i„'s are the homogeneous symmetric functions; /i„ = for n < 0. By the 
Cauchy-Binet formula the last sum is readily seen to be equal to 



det 



hk-\i+i{h)hk-vj+j{(' 

lk>-l 



and the sum inside the determinant is equal to the needed matrix element of L 

because for p,q > —I we have 

X] hk-p{b)hk-q{a) = ^hk-p{b)hk-q{a) = ^hk+q-p{h)hk{a) 



k>-l 



and 



^^hk+n{b)hk{a)C = E f'kAb)& E f'kMC' 



expE(p„(a)C-"+p„(6)C). □ 



It is convenient to introduce separate notations for the two values of the kernel 
entering the statement of Theorem 2.2: Define Kj^[a,x;T,y) and K-{a,x;T,y) 
through the generating functions 

We will denote by K± [a, r] the Z' x Z' matrices with matrix elements 

{K±[a, T]):^y = K±{a, x; r, y). 

Here a and r are allowed to take values between and N. 

Lemma 2.10. (i) For any a G {0, . . . , A''} we have [a, a] = 1 + K_ [a, a] . 
(ii) For any a,T & {Q, . . . ,N — 1} we have L[a + 1]K±[<j,t\ = K±[a + 1, r]. 
(Hi) For any t G {0,. . . ,N} we have QK_[N,t] = -K+[0,t]. 

Proof. All these statements are proved by simple algebraic manipulations. 
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For (i) we have 



meZ' 



m£Z' 



whence K+[(j,(j] — K-[(j,a\ = 1. 

The formula (ii) follows from the definitions of i[cr], K±[a, r], and the fact that 

(exp E + 1]C-" + bn[cJ + 1]C))F{^, C) = + 1, 0- 
Finally, for (iii) we have 

zF{N,tC) ^ (iCr?)'" 



i:«3A-_|iv,x|,..CV = -f^Ei 

F 



The relations of Lemma 2.10 immediately imply that if we introduce a X x X 
matrix K which has the block form 



K = 



i^+[0,0] K+[0,1] K+[0,2] 

K_[1,0] K+[l,l] K+[l,2] 

kJ[2,0] kJ[2,1] K+[2,2] 

.K_[N,0] K_[N,1] K_[N,2] 



K+[0,N] 
K+[1,N] 
K+[2,N] 

K+[N,N]. 



(the (cr, T)-block is equal to i('+[iT, r] if a < r, and to K-[<j,t\ otherwise), then 
we have the matrix relation (1 + L)K = L. In order to extract the probabilistic 
meaning of this relation we need to introduce certain finite point approximations 
of the shift-mixed periodic Schur process. 

For any m = 1,2,... denote by U{m) the subset of Z' consisting of 2m half- 
integers situated symmetrically around 0: 



S'(m) = {-m-l- 5,-m + 



3 _i 1 

2'"""' 2'2'""'' 



m- |,m- i}, 



Denote also by L(™)[cr], K^^\a,T\, and Q^"' the restrictions of the Z' x II 
matrices L[cr], /^±[cr,r], and Q to Z'(m) x Z'(m), and denote by L^™) and 
the block matrices built from L^^^fa], ^ST^"^ [tr, r] , Q^'") in the same way as L and 
K are built from L[cr], Jr±[cr, r], Q. Set 



X^") = Z' (m) U • • • U Z' (m) ( + 1 copies) . 



Note that 



3e(™)xje('") ' 
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Lemma 2.11. (i) det(l + L^™)) = z'^t-'^{l + 0(1)) where 0{1) stands for an 

element ofAhoii^)- 

(a) All matrix elements of (1 + L'^™')^^ belong to Ahoi- Consequently, for any 
X c we have deti^Vdet(l + L^™)) e Ahoi- 

Proof, (i) Since det(l + i^™)) = Y^xcx'-^'> det-£/l^\ we need to see which minors 
det yield terms of the lowest possible degree. Since 

{L[k])^y & Ahoi{\x - y\) and (i[fc])^^ = 1 + 0(1), k = l,...,N, 



it is immediate that the only lowest degree term comes from det i^"* with 
X= {-m+i,...,-i}u-..u{-m+i,...,-i} 

(TV + 1 copies of the same set of negative elements in Z'(m)), and it is equal to 

z"^t~~ (recall that z does not contribute to the degree count). 

(ii) Matrix elements of (1 + i^™))"^ are ratios of the linear combinations of 

minors of L^™) and dct(l + i'™)). Since all minors of i^™^ lie in t~^ ■ Ahoi, by 
(i) we see that the matrix elements belong to Ahoi- 

As for the ratios det Lx / det(l + -L^™)), observe that they are linear combinations 
of ratios of the form det(l + L("*^)x/det(l + L^'")) which coincide, up to a sign, 
with minors of the inverse matrix (1 + L^™')^^. □ 

We arc now in a position to prove that L^^^-enscmblcs on X^™-* approximate the 
shift-mixed periodic Schur process on X as m becomes large. 

Proposition 2.12. The values of the correlation functions of the L^"^^ -ensembles 
restricted to the last N copies ofl'i^m) in X(™) =Z'(to)U---UZ'(to) converge, as 
m — > DO, to those of the dynamical correlation functions of the shift-mixed periodic 
Schur process in the topology of A. 

Proof. In order to prove this statement we will construct an injective map that 
associates to any point configuration X C X'-™^ of nonzero weight in the i^™^- 
ensemble a sequence of partitions {\^^\ . . . , A*^^)) and an integer S so that: 

• The intersection of the point configuration 

{5 + A«-^ + i}^>^U...u{5 + Ar)-.+ i}^^^ 

with Z'(m) U • • • U Z'(rn) coincides with the restriction of X to the last N copies of 
Z'(m) in X^™). 

• The weight of X in the i^'^^-ensemble and the weight of (A, 5*) in the shift- 
mixed periodic Schur process are obtained from each other by multiplication by a 
constant of the form Cm = 1 + 0(1) and by addition of an element of a high enough 
degree: 



detLl^) 
det(l + L^^) 



Cm ■ z^t^ J2 m) + 0{d{m)), 



where d{m) does not depend on X and d{m) — > oo as m — > oo 



10 



^"We should have used Cm/{9:i{z',t)Z{N,t,a,b)) instead of Cm thus taking into account the 
partition function. However, this is equivalent because 03{z; t)Z{N, t, a, b) and l/{03{z; t)Z{N, t, a, b)) 
are both of the form 1 + 0(1). 



PERIODIC SCHUR PROCESS AND CYLINDRIC PARTITIONS 



23 



• The lowest degree of the weights of pairs (A, S) not covered by this map goes 
to infinity as m — > oo. 

Since ah the weights add up to 1 in both the L^^^-ensemble and the shift-mixed 
periodic Schur process, the existence of such a map imphes that the degree of Cm — 1 
goes to infinity as m — > oo, and the needed convergence of the correlation functions 
readily follows. 

Observe that if a set X C X^™-* has a nonzero weight in the L^^^-ensemble (i.e., 
det l'"^^ 7^ 0) then its intersections with all + 1 copies of Z'(m) in X'-™-' must 
have the same cardinality because of the specific block structure of L^"^\ and its 
intersections with the first and the last copies of U{m) must coincide because the 
matrix Q is diagonal. Thus, without loss of generality we may assume that 

X = . . .,xf^) u • • .,4'} u • • ■,xf'} u • • • u . . 

for some I < 2m, where xf'^ = x\^^ for all 1 < i < Z. Our notation for X means 
that the first group of points lies in the first copy of Z'(m), the second one lies in 
the second copy of Z'(to) and so on. 

The element (A^^^ . . . , A^^^; 5*) of the shift-mixed process corresponding to such 
X is defined as follows: S = I — m and 

{5 + Af)-i+i}.>,={4^\...,4^)}u{-m-i,-m-i,-m-|,...} 

for any t = 1, . . . ,N. It is readily seen that this formula correctly defines the 
partitions A^^). Note that ^(A^^)) < I for all r. 

The pair (A^^^ . . . , A^^^; S) is not covered by this map if and only if for some 
r the point configuration {S + x'f '' — i + ^}^^-^ either does not contain the set 
Z'^_^ or has a nonzero intersection with Z>^. In the latter case we must have 

Xi^^ + S > m + I and in the former case we must have {X^'^^Yi — S > m + 1. This 
means that either l^l or lA^'^)] is > 

The definition of the shift-mixed process implies that the weight of (A, S) has 
degree at least 

f + ||A«| - |A(2)|| + .. . + ||A(^-i)| - |aW|| + |aW| < f + ^max^|A(^)|. 

Therefore, the pairs (A. S) not in the image of our map have weights of degrees 
uniformly going to infinity as m ^ oo. 

It remains to compare the weights of a point configuration X C X*-™^ in the 
X('")-ensemble and of its image in the shift-mixed process. 

Using Lemma 2.11(i) we obtain that the weig ht detL^Vdet(l + L^")) of X 
equals 

const -^-"t^ det[(L[l])^a)_,(o)];_ • • • det[{L[N])^^.,^^(.-^,][.^^ [] zf^^"'. 

i—l 

with const = 1 + 0(1). Collecting the powers of z and t yields 
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where we used the definition of S = l-m and the fact that ^(A^^)) < Evaluating 
the minors of L[t\ by Lemma 2.9 (which is applicable because < I) we obtain 

det[(L[l])^(i) (0)]' • • • det[(L[7V])^(«) = const 

^ s;,(0)/^(i) (a[l]) (&[!]) + 0{2l - £(A(°)) - ^(A^i)) + 2) 




J2 Sa(i)/^(^) (a[2]) (6[2]) + 0(2Z - i{X^'^) - ({X^'^) + 2) 



Yl «A(--)^(-) («[^]) Sa(")/^(") m]) + 0{2l- £(A(^-i)) - ^(AW) + 2) 

with const = 1 + 0(1). To conclude the proof we need to show that all the re- 
mainders O ( • ) in the expression above can be removed by adding the correction of 
degree uniformly going to infinity as m oc. 

Observe that the degree of the factor 'I is bounded if and only if \S\ = 

\l — m| and |A^^^| are bounded. Then in the last factor in the product above the 
term 0(21 — i{X^^~^^) — £{X^-^^) + 2) has degree going to infinity with m (because 
I - £(A(^)) >m-S- I A(^) I), so if we want the weight of X or the weight of (A, 5) 
to be of bounded degree, the degree of the sum over /x^^^ has to be bounded. Since 

degY^sx/^{a)s^/i^{b) > ||A| - 

this means that jA^-^"^)] has to be bounded. Repeating the argument with the 
second to last factor we conclude that |A(^~^^| is bounded, and so on. 

The final conclusion is that for the corresponding X and (A, S) such that the 
minimum of the degrees of their weights is bounded, we must have |Z — m| and 
lA^"^)], r = 1,...,A'', bounded. But then all C(-)'s in the product above have 
degrees uniformly going to infinity as m ^ oo, and thus the difference of the 
weights of X and (A, S) has degree uniformly going to infinity as m ^ oo. The 
proof of Proposition 2.12 is complete. □ 

Let us now conclude the proof of Theorem 2.2. From the determinantal formula 

for the correlation functions of general L-cnscmbles, sec the beginning of the proof, 
we know that the correlation functions of the L^™) -ensembles are given by minors 
of the matrices L^"*' (1 + L^"^^)~^. The last step of the proof is to show that matrix 
elements of these matrices converge to those of the kernel K in the topology of A. 
The definitions of matrices L[a] and K[(7,t] imply that 

deg{L[a])ooy>\x-y\, iK[a,T])^y > \x - y\, x,yeZ', 

for any a and r. Hence, from Lemma 2.10(ii) we obtain {x,y € Z'(m)) 

(L(™)[a+l]ifi")Kr])^^ = {K^^\a + 1,T]U 

+ C(min{m + ^ — x,x — m — ^} + min{m + ^ — y, y — m — ^}) . 



^This is the computation mentioned in Remark 2.7 above. 
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The relations (i) and (iii) clearly remain unchanged when restricted to Z'(m): 
Hence, 

((1 + = 4™) + 0(min{m + 1 - y, y - m - ^}) . 

Multiplying by (1 + L^"*^)"^ on the left and using Lemma 2.11(ii) we obtain 

= ((1 + + 0(min{m + ^ - y, y - m - ^}). 

Since min{TO + ^ — y, y — m— ^}— >ooasm— >oofor any fixed y, the needed 
convergence follows. The proof of Theorem 2.2 is complete. □ 

3. Bulk scaling limit 

As was mentioned in §1, the presence of the parameter t is crucial for defining 
the periodic Schur process; at t = 1 the partition function is, generally speaking, 
infinite. Thus, one might expect that if t ^ 1— then the random Young diagrams 
become large. Indeed, this is correct: In the analytic setting we will show that as 
t 1—, the density fimction of the scaled random point configuration | lni| ■C{X^'^^ ) 
tends to a nontrivial limit, and this limit is independent of r. The main result of 
this section is the computation of the local limit of the correlation functions of the 
periodic Schur process and its shift-mixed modification near points of fixed global 
limit density. 

Denote r = \nt~^. Throughout this section we assume that r > (equivalently, 
< t < 1), and also that the convergence conditions of Remark 2.3 are satisfied. 
Namely, we assume that |a„[A;]|, \ bn[l]\ < const -R" for some < i? < 1 and all 1 < 
k,l < N, and we also assume that | arg(^;)| < tt. 

Recall the notation A). = a/s[?Ti], B). = J2m=i ^k[m], introduced in §1. 

Theorem 3.1. (i) Assume that Ak = Bk for all k = 1,2, ... . Then, as t ^ 1, 
the dynamical correlation functions of the shift-mixed periodic Schur process have 
a limit in the following sense: Choose xi{t), . . . ,Xn{t) € Z' such that as t ^ 1, 

rxk{t) 7 for all k — I, . . . ,n and some 7 G M, and, all pairwise distances Xi — Xj = 
Xi{t) — Xj{t) are independent oft. Then for any 1 < n, . . . , r„ < iV 

limpf ft(Ti,ari(i); . . . ,T„,a;„(i)) = det [4^;;!:) (x, - xj)]; .^^, 
where the correlation kernel has the following form 



ICi^''>\d) 



|^£l| 1+^-1 exp(7-i:„>i(A„f-"+B„C")) ^ - ^' 
1_ r exp(i:„> i:g^^^,(an[fc]C-"+bn[fc]C")) dC 



(ii) Under the same assumptions the dynamical correlation functions of the periodic 
Schur process converge to the limiting expression above evaluated at z = 1: 

lmipn{Ti,Xi{t); . . .,T„,Xn{t)) = dct [4;?] (Xj - Xj)]"!.^^. 
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Comments. 1. The limit correlation functions as functions on ZU • • • UZ {N copies) 
are invariant with respect to the simultaneous shifts of all variables. 

2. It will be clear from the proof that a slightly more general statement is true: 
For several groups of variables {x^k^\t)}'^^-^, m = 1, . . . , M, such that rx^^^ jm, 
the pairwise distances inside every group are independent of t, and the distances 
between different groups tend to infinity: 

the limit of the dynamical correlation functions of the shift-mixed process is the 
product of determinants: 



M 



limp^ ({r(™\a;l'")(i)}i<,<„„,i<„<M) = J] dei[K^:dTlU^r ' 

m=l « ' J 

and the same is true for the correlation functions of the periodic Schur process with 
= 1 in the right-hand side of the formula above. Roughly speaking, this means 
that particles become independent as the distance between them grows. 

3. The global limit density function mentioned in the beginning of the section is 
equal to 



27riJf^l=il + ^-ieT-S„>x(^nC-"+B„C'') ( 

for the shift-mixed process, and to p*-^-* (7) for the unmixed process. This formula has 
the following corollary: If one assumes the existence of the limit shape, as t ^ 1, 
of the random Young diagrams A^"^^ distributed according to the periodic Schur 
process, then this limit shape can be easily guessed. Denote by i and j the row and 
column coordinates on the diagrams and introduce new coordinates u = r{j — i) 
and V = r{i + j). Then the equation for the boundary of the hypothetical limit 
shape has the form 

v{u) =u + 2<f In f 1 + e-«+Sn>i(^nC-"+B„C'')\ ^ . 

7|ci=i ^ ^ C 

This formula follows from the relation p{u) = ^(l — v'{u)), see [BOO, Remark 1.7] 
for an explanation. 

4. In the case of the uniform measure on partitions, when the specializations 
a[fc],6[A;] are trivial. Theorem 3.1 (or rather its extended version from Comment 1 
above) coincides with Theorem 7 of [01]. In this case the limit correlation kernel 
degenerates: 

6od-{l + z'^e-<)-', a<T, 

-5od • (1 + ^e-T)-\ a>T. 

For the uniform measure the limit shape does exist, see [V], and the formula for 
v{u) from the previous comment produces the correct answer: v{u) = u + 21n(l -|- 

e~") or + = 1. 

Proof of Theorem 3.1. (i) We start with the integral representation of the correla- 
tion kernel for the shift-mixed process from Remark 2.4. Replacing the integration 
variable by ^ = C»? we obtain 

K( \- nn>i(l-^")' / / F{a,0 03iz^;t) dCd^ 

K[a,x,T,y) g^^,.^t){2mr J J F{T,a-') es{-^t-i;t) ^-v+^^v+i ' 

|C| = l|a=l±e 
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where e > is much smancr then r. and we choose |^| = 1 + f for a < t and 
1^1 = 1 — e for cr>T. The next step is to fix an arbitrary C, on the unit circle and 
to evaluate the asymptotics of the integral over ^. 

Proposition 3.2. (i) Assume that \ &vgz\ < n. Then, as t ^ 1, the function 

n„>i(i-i")'^3(^?;t) 



/(0 = -- 



on the circle |^| = lie uniformly converges to on the complement to any neigh- 
borhood of the point ^ = 1 ± e. On the other hand, there exists 6 > such that for 
I arg(C)| < S 

HO = .Ja ^ • (1 + /o(0) 

' e r- — e 

where /o(C) an analytic function which, ast^ 1, uniformly converges to while 
''/o(^) remains uniformly hounded}"^ All the estimates are uniform, in z varying in 
a compact set of the complex plane bounded away from the negative real semiaxis. 

(ii) Assume that z lies on the unit circle \z\ = 1. Then on the circle |^| = 1 ± e, 
as f — > 1, we have the bound \ f{C)\ ^ const -e"^ which is uniform in both z and ^. 

Proof, (i) Denote v = ^^^^ ^'^d w = 2^1n2:. Applying the imaginary Jacobi 
transform, see e.g. [Er, 13.22(8)], we obtain 

6i{z;t) = {^Y e —6^{e—;e- — ), 

e^{zt t) = (f ) ^ e '-^03 (e^^ ; e"^) , 



Hence, 



We have 



^ n>l ^ 



^3 e r ; e r (^3 -e r ; e 



iini!2^rr(i-rf = i. 



The last relation can be obtained, for example, from the imaginary Jacobi transform 
oi di{x;t) because 9i{x;t) ~ {x^ — x~^) {t;t)l^ as a; — > 1. Thus, 



27r^u(l-2m) 



_ (^) ^ J] (1 _ t^f e-^((''+-)^-^-(-iJ+i)=) = M£ : ""(i + ^W) 



n>l 



^Here we use the principal branch of the logarithm function for both ln(z) and In(^). 
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Writing down the products for ^3's explicitly we obtain 

93(6 r -e r j =Uo- [[ il + e r 1(1+6 " ) 

^3(e^^;e- — J =no- M + e fl + e ^ 

n 1 + e ^ ^ 1 + e F 



X 

n>0 



^sf— e r ;e '^l=no-(l — e*- 1 

n( 47r2(„-n-l) \ / 4,r2(-t,-n-l) 

II — e »■ 1(1 — e »■ 



X 
n>0 



where Ho = nra>i(l ~ ^ ^ 1 as r ^ 0. Using the principal branch of the 

logarithm we may assume that 

= ^arg^ e and Jftw = ^ argz e (-i, i) 

(recall that | argz| < tt by assumption). 

It is readily seen that under these conditions, as r ^ 0, all three products over 
n > in the formulas above uniformly tend to 1 while their derivatives with respect 
to V multiplied by r are uniformly bounded. The remaining factors are 



27ri e ^ ( 1 + e ? — - ) ( 1 + e" 

/(O 



4„2(„_^„_l)\ / 4,r2(_„_„_l) 



Assume first that is bounded away from (that is, ^ is bounded away from 
^ = 1 ± e). Note that since e <C r, we have SJt; = ^ In |1 ± e| = o(r) and 

SR(u(l - 2w;)) = m.v(\ - mw) + o(r). 

If < and 5i(u -\-w) > — ^ then the absolute value of the asymptotic expression 
for /(^) above is bounded by 

const •r-ie^-'^-d-^^-)'^"' 

which uniformly converges to as r — > 0. li < and ^{v + w) < —\ then the 
bound is 

const • r-^e^''^^^^''^^"^*^"''"^*^''"^*^'""^^'''' = const • r-ie~^''^(*^''+^^^^*^'"+^^'~~' 
which also goes to as r ^ 0. If 3?t; > and "^{v + w) < ^ then |/(^)| is bounded 

by 

const •r-ie2-'(*^''(i-2*^-)-2*^'')'-"' = const •r-ie-2-'*^''(i+2*^«')'--' 
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which goes to as r — > 0. Finally, if^v > and + w) > ^ then the bound has 
the form 

const .r.-ie2''^(*^''(^~^*^'")+^(*^''+*^'"~5)~^*^''^''~' = const .^-le"^''^'*^''"^)'^"^*^'"^''"' 

which is also small as r — > 0. This takes care of the first statement (about uniform 
convergence to 0) of Proposition 3.2(i). 

Since |5ftw| < i, we can choose S > such that |3??/) i ^1 < 5- Then for ^ 
with I arg(^)| = 27r|3?w| < S we have ^{v + w — ^) < 0, and the two factors in the 
numerator of the approximation for /(^) above uniformly converge to 1 as r ^ 
with their derivatives with respect to v multiplied by r are uniformly bounded. 
Thus, 

/(?) ~ / -2.2., 



r e 



as desired. It is also immediately visible that all the estimates above are uniform 
in z varying in a compact set bounded away from the negative semiaxis. The proof 

of (i) is complete. 

The proof of (ii) follows the arguments used above to prove the first part of (i). 
There are two differences in estimates: the range of is now the whole segment 
[— ^,^] rather than a closed subset of the interval (—5,5), and since is not 
bounded away from 0, we have an additional potentially small factor of the form 
in the denominator. We have 

arg(e±^) = ±^ = ±^ ln(l ± e) ~ ±^ 



and hence 



1 - e 



> const Adding this estimate to those derived above 
yields the needed bound. □ 

Our next step is to compute the asymptotics of the ratio F{a, () /F{t, C4^^) from 
the integral representation of the correlation kernel given in the beginning of the 
proof of Theorem 3.1. Recall that the definition of the function F{t, was given 
just before Theorem 2.2. After simple manipulations, using the assumption that 
B/j = Afc for fc > 1, we obtain 

Ca N T N \ 

fe=l fe=cr+l k=l k=T+l ) 

For any a and t this expression viewed as a function in (C,C) ranging over the 
circles \C,\ = 1, |^| = 1 ± e, remains uniformly bounded away from and 00 as f — > 1 
if < e ^ r. Indeed, the boundedness of the second factor is obvious, while for the 
first factor we have 

1 - r > nrr, 1^=^" - = 1(1 ± e)=^" - 1| < ne{\ + e)" 
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and hence 



E 

n>l 



3? 



A,;"(i-r")--ur"(i-4") 



E 

n>l 



3? 



1 -t" 

< 2er-^^i-"(l + e)"|A, 



n>l 

is bounded as t — > 1 because er~^ — > 0, and \An\ = 0{R") for a fixed -R < 1 as 
n ^ oo. Similar arguments imply that the derivative ^ p^^^'^^i^ multiplied by r is 
uniformly bounded. 

Let us now look more carefully at the ratio in question when ^ is close to the 
real line. Assume that | arg^| < a with r <C a <C . Then, using the estimates 

|^±" - lipnln^l < (n| ln^|)2max{l, 1^1"} < const •n2a2(l + e)"«:nV(l + e)", 

1 



(1-r) 



n\-l 



nr 



It" - 1 - nr\ ^ n^r^ 



(1 — t") nr nrt^ nr 
and the fact that |j4„| = 0{R"') as n — > oo, we obtain 



= t- 



= (1 + ^o(O) • exp(r-i In^ + ^nC"")) 

-t^lT^'S,? J n>l 

/a T N 

X exp^ r"E""W-^""E""[^]-^" E ^»w+c" E ^"t'^i 



AT 



n>l 



fe=l 



fe=l 



fe=o-+l 



fe=T + l 



where -Fo(0 uniformly converges to as r — > 0. Differentiating this relation with 

d€ ^'(t,«-i) 



respect to ^ and using the boundedness of ^(^^^^-i) and r^^ pf}rf-i\ well as the 



fact that 



n>l 



exp(r-Mn$^(^„C" + ^nr")) = exp(r-Mn(l±e) + A„r")) ^ 



n>l 



as r ^ 0, we see that r-FQ(^) is uniformly bounded in both ^ ranging over the arch 
I arg^l < a, 1^1 = 1 ± e, and C ranging over the unit circle. 

We can now proceed to evaluating the asymptotics of the integral over ^. Set 



n\3 



G±(C) = - 



F{a,0 e3{z^;t) 



|{|=i±. 



According to the hypothesis of Theorem 3.1, we assume that y = 'yr~^{l + o(l)). 
Note that l^*'"'' = | ^ 1 as r ^ 0, because er 0. 

Let us split the circle |^| = 1 ± e into three parts: | arg^| > S with ^ > taken 
from Proposition 3.2(i); a < \ arg^j < S with r <C a ^ ; and | arg^| < a. 

The integral over the first part tends to zero (uniformly in on the unit circle) 
because /(^) of Proposition 3.1 converges to and the factors (^^^^-i) and 
remain bounded. 
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The integral over the second part, by virtue of the boundedness of the same 
factors and the asymptotics of /(^) from Proposition 3.1(i), is bounded by (recah 
the notation v = ^^^S.) 



const 



e 



e r 



dv. 



Introducing a new variable v = 27r^r ^v, we can rewrite this expression as 
const / 



e ^ 



dv, 



which tends to as r — > as long as | argz| < n, er ^ 0, and a/r — » oo. 

The only nonzero contribution comes from the third part. Using the variable 
ij = 2n'^r^^v = —iT:r~^\nf. the estimates for the ratio „f^'^'S}i\ obtained above, 
and Proposition 3.1(i), we can write the integral over the third part | arg^| < a in 
the form (note that = ^dv) 

Co- r N N \ 

fe=l fe=l fe=cr+l fe=T+l / 



1 r „f (i"W-7+En>i(>inr+Anr")) 

J {l + G,iv))dv, 



X — 

ITT 

where the function Go{v) uniformly converges to as t ^ 1, and the derivative 
Gq{v) remains uniformly bounded. The statement of Theorem 3.1(i) is a corollary 
of this formula, the above estimates, and the following lemma. 

Lemma 3.3. For any a € C such that |3a| < 1 the following limit relation holds: 

f +oo±i£ 



/ 7i — 7^ + x))dx 



1 + e^'^°- 



where we assume that the function g{s, x) on M ± is uniformly tends to as e ^ 

0+, while its derivative -j^g{e,x) remains uniformly bounded. The convergence is 
uniform in a varying over any compact subset o/C with max|3a| < 1. 

Proof. The results with x &M. + ie and x &M. — ie are obtained from each other by 
the change of sign of the integration variable. Hence, we may assume that x varies 
over R — ie. 

Let us handle the term with g{s, x) first. Split the integral into three parts: 

/■°°-'^ e'''-='g{x,s)dx _ r-^o-ie ^ jx„-re ^ .+<x>-fe 
J — oo—ie ^ ^ J — oo—ie J —XQ—ie J XQ—ie 
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for some a;o = a:o(e) G R that will be chosen later. Denote M{e) = sup^^gj^.j^ \g{x, s) 
The third integral can be estimated as follows: 



+00 — ie iax 



e^°'^g{x, e)dx 



— e 



< 



M{e) f 

J X, 



M{e) f 



+00 ^-Xa)vdy 

ev - e-y 



+C30 — ^e 



M{e) 



1 - 9a 



In 



g(l-|9!o|)a;o 
g(l-|3;o|)a:o _ 



If Xo(e) ^ as £ 0, then the last expression is bounded by const ■M{e) \ lna;o(£)|, 
which tends to if Xi^{e) > e~^^ . A similar estimate holds for the first integral. As 
for the second integral, using the notation M'{e) = sup \-^{^^°'^9{x,£))\, 



{—xo — ie,xo — ie) 



we have 



XQ — ze 



oX — X 



/xo 
-Xo 



\y\dy 



Igy-fe _ g-y+iel 



g{-is,e)dy 



< const M'{e)xo 



ln(e^ 



1) - ln(e 



1) - In ■ 



_ g-y+is 



< const (M'(£)a;o + M{e)) 



which tends to if x{e) — > as £ — > 0. Thus, we proved that the term with g{x, e) 
tends to as £ ^ 0. 

The remaining term gives (we use y = cosh(?/) / sin e below) 



e-rfx r 



dy 



gy-is _ g-y+is 



+ ie COS £ 



cos(ay)(e^ + e~y)dy 
{ey - 6-2^)2 + 4 sin^£ 

cos(a • arccosh(y sin£))dy 



2e^0j_ 



sin(ay)(e^ — e y)dy 
{ey - e-yy + 4sin^£ 
sm{ay)dy 



ey - e-y 



ITT 



ITV ITT e" 

Y + Y + 1 ~ YTe-^ 



a 



Lemma 3.3 implies that as t — > 1, 



AT 



G±(o-±exp^ r"E"»w-c""E«"W-c" E bn[k]+c fe^w 

rt>l V k=l fe=l fe=cr+l 

1 

X 



N 

^ E 

/s=r+l 



1 + g=F(lnW-7+E„>i(^nC"+A„f-")) 



where the convergence is uniform in ^ varying over the unit circle, and in z varying 
in a compact subset of C not touching the negative real semiaxis. Since 



K{a,x;T,y) 
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this completes the proof of part (i) of Theorem 3.1. 

In order to prove Theorem 3.1(ii) we will use Proposition 2.1, which (along with 
Theorem 3.2) implies 

Pn{ri,Xi;...;Tn,Xn) = 77^/ <ie^[K{Ti,Xi;Tj,Xj)]ij^i93{z;t) — . 

Using the familiar notation w = ^Inz = ^argz G [— ^;^] and applying the 
imaginary Jacobi transform, we obtain 

(this formula was already used in the proof of Proposition 3.2 above). The product 
formula for ^3(6^^^ ;e^^) implies that this function of w remains uniformly 
bounded on [— |, |] as r ^ 0, and it uniformly converges to 1 on any open interval 

inside [—^ ^ ^ 



2' 2J 



Let us split the domain of integration over z in the formula for /9„ (n , a;i ; . . . ; r„ , a;„) 
above into two arches: | arg2;| < c < tt and | arg2;| > c, where c > in a arbitrary 

constant < tt. 

The integral over the first arch, by Theorem 3.1(i), is equal to 



1 /■' 



detitcCr-y^Xi - xj)]l.^^ dw + 0(1) 



as r — > 0. Since (^) ^ e~ ^"r" converges to the delta-function at t« = as r ^ 0, 
the integral above converges to det[ici]'^] {xi — Xj)]^_^_^. 

As for the integral over the arch | argz| > c, we use Proposition 3.2(ii) and the 
boundedness of the ratio F{a, Q)/F{t, C^~^) proved earlier to see that the absolute 
value of the integrand is bomidod by const -e r 2r, which converges to as 
r ^ as long as In <C r. Since this does not contradict our previous assumption 
that e <C r, we may ignore the integral over the second arch in the limit r ^ 0. 
This completes the proof of Theorem 3.1. □ 

Example 3.4. Consider the periodic Schur process with N = 1 and 



ln[l] = bn[l] 



^, n = l, 
0, n > 1, 



with an arbitrary > 0. This is equivalent to considering a probability measure 
on pairs (A D n) of partitions given by 

ur^ 1 ( TT 1 \ 'dim2(A/M)cl^lcl''l 

Prob{A,At} = ei-<=i<=2 — ^ , 

V i>i -(^1^2)"/ (1^1 -ImI)!' 

where Ci = t'&'^, = and dim(A//u) is the number of standard Young tableaux 
of shape A//i. Here we used the relation 

( dim(A//x)^l^l-l^l 
sx/M^\) = sx/MA) = — (|;^|_|^|), — 
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and the formula for the partition function from Proposition 1.1. 

Theorem 3.1(ii) yields the following limit result for the local correlation functions 
of this measure: If a;i(t), . . . , .T„(t) <E Z' are such that 

lini (ini^^ • a;i(t)) = 7, i=l,...,n, 

and the pairwise distances Xi — xj do not depend on t, then 

liniProb{{a;i(t), . . . ,a;„(i)} c{Xi-i+ i}i>i} = det[IC^^\xi - 

where 

1 dC 



IC=i| 

The formula for the limit correlation functions for /i is exactly the same. 

In the limit •& the measure on A's approaches the uniform measure Prob{A} = 
const -fl^l on partitions, and, correspondingly, )C^'^^{d) — > 5od • (1 + cf. Com- 

ment 4 to Theorem 3.1. 

As 00, the measure on A's looks more like the poissonized Plancherel mea- 
sure on partitions with weights of the form Prob{A} ~ dim^ A • constl'^l/|A|!^, be- 
cause as C2 — 'd~'^ 0, the partition /i tends to be small. When 7 is of the same 
order as the local correlation kernel converges to the discrete sine kernel, 

r 0, e > 2, 

r y-(-y)rj\ ^ I ^2\uc\>e dC I sin(arccos(6/2) • d) 

Sod, e < -2, 



i9— >oo 

i/^^e |c=i| 



which is exactly the bulk scaling limit of the correlation kernel of the poissonized 
Plancherel measure with large poissonization parameter, cf. [BOO, Theorem 3]. 

4. Extensions of the discrete sine kernel 

The goal of this section is to construct an infinite-dimensional family of determi- 
nantal point processes on Z x Z such that the restrictions of their correlation kernels 
to copies of Z obtained by fixing the first coordinate coincide with the discrete sine 
kernel. 

Two such extensions of the discrete sine kernel have been constructed previously: 
in [ORl] a kernel called the incomplete beta kernel was obtained in the bulk limit 
of the large imiformly distributed plane partitions, and in [B02] another extension 
was obtained in the bulk limit of the Markov chains on partitions preserving the 
Plancherel measure (equivalently, in the bulk of the multi-layer polynuclear growth 
process with droplet initial conditions). Both these extensions are included in the 
family that will be constructed below. 

We start with some necessary generalities. 

Set X = Z X Z and let 2^ be the set of all subsets of X. The first coordinate of 
points in X will be viewed as "time" while the second coordinate will be viewed as 
"space" . Pick any exhaustion of X by a sequence of increasing finite sets: 

BiCB2CB3C---cX, |-B„|<oo, \jBn = X. 

n>l 
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For instance, one can take BnS to be growing boxes B„ = {{x, t) : \t\ < n, \x\ < n}. 
Clearly, 2^ = Km 2^", and we equip 2^ with the projective limit topology. In other 

words, a sequence C 2^ converges to X C X if and only if for any m> 1 there 

exists M = M{m) such that for n > M we have Xn n Bm = Xf\ Bm- It is easy to 
sec that the topology does not depend on the choice of BnS. Since projective limits 
of compact topological spaces are compact, the space 2^ is compact with respect 
to this topology. 

A random point process on X is, by definition, a Borcl probability measure on 
2^ 13 Qjjg construct random point processes on X is to provide a count- 

able set Xi,X2,... of subsets of X and the set of positive weights wi,W2,---; 
ProbjXi} = tUi, satisfying the normalization condition X^j>x Wi = 1. 

The correlation functions p„ of a random point process on X are probabilities 
for random subsets of X to contain a given finite set: 

pn{ti,x-i; ■ ..;tn,Xn) = Prob{X C X | {(ti,xi), . . . , (T„,a:„)} C X}. 

For a finite set B C X, the values of the correlation functions on B define the 
projection of the measure on 2^ to 2® uniquely. Indeed, the two sets of 2l^l — 1 
numbers giving for any nontrivial subset Bq of B the probability that the intersection 
of the random set X c X with B contains Bo or coincides with Bo are related 
by a nondegenerate linear transformation obtained from the inclusion-exclusion 
principle. One set consists of the values of the correlation functions while the other 
one consists of the weights of the subsets with respect to the projected measure on 
2^. For example, for B = {bi, 62} we have 

Prob{X e 2^ I X n B = {61}} = - P2(6i, 62), 

Prob{X e 2^ I X n B = {62}} = Pi (62) - P2(6i, 62), 
Prob{X e 2^ I X n B = {61, 62}} = P2{bi, 62). 

Lemma 4.1. Let Vi,V2, ■ ■ ■ be a sequence of random point processes on 1? , and 
assume that, as m —> 00, all correlation functions of Pm converge pointwise: For 
any n = 1, 2, . . . and any {Ti.Xi) G I? , i = 1, . . . , n, we have 

lim pn{Ti,xi] ...]Tn,Xn \ Vm) =■ r„(Ti,xi; . . . ;t„,x„). 
Then there exists a unique random point process V on 1? such that 

Proof. The uniqueness of V was demonstrated above. The same argument shows 
that for any finite set B, the set of possible values of the correlation functions of 
a probability measure on 2^ is described by a finite list of linear inequalities. For 
instance, in the case B = {61,62} these inequalities are 

Prob{0} = 1 - pi (61) - pi (62) + P2(6i, 62) > 0, 
Prob{6i} =pi(6i) -p2(6i,62) > 0, Prob{62} = ^1(62) - P2(6i, 62) > 0, 
Prob{6i,62} = P2(6i,62) > 0. 



^^In the conventional terminology, this is actually the definition of a simple or multiplicity free 
random point process. Since all our processes axe multiplicity free, the definition above is general 
enough for our purposes. 
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Clearly, these inequalities are preserved under limit transitions. Hence, for any 
of the finite sets Bm used above there exists a unique probability measure Pm on 
2^"* such that its correlation functions coincide with the restrictions of the limit 
functions r„ to Bm- Moreover, these measures are consistent: for mi > TO2 the 
projection of the measure Vmi on 2^"i to 2^"'2 coincides with ■ (Indeed, both 
these measures have the same correlation functions r„|B„^.) Since 2^ = lim2^'", 
we can take V = limP^. □ 

Observe that all that was said above is not specific to the set X — "Z'^ and applies 
equally well to any discrete countable set X. In particular, we will consider random 
point processes on Xn '■= {1, • • • , N} x Z below. Since Xn C X, a random point 
process on Xat may also be viewed as a process on X whose point configurations 
are contained in Xjv almost surely. 

Let us now return to the periodic Schur process. So far we have not discussed 
the positivity of the weights which we used to define the process. One way (but 
not the only way, see e.g. §8 below) to guarantee the nonncgativity of the weights 
W{X,fj,) introduced in the beginning of §1 is to take t > and to demand that 
all the specializations of the skew Schur functions are nonnegative. The following 
classical result is useful. 

Due to the Jacobi-Trudi formula 



s\/u. = det[hx,-i-^,^+j]lj^i, r > max{£(A),£(/x)}, 

see [Macd, §1 (5.4)], we need to guarantee the nonnegativity of the determinants in 
the right-hand side (here hnS are the complete homogeneous symmetric functions). 
Recall that a sequence {c„}^q is called totally positive if all minors of the matrix 
[ci-j]ij>o are nonnegative. Here all C-k for k > are assumed to be equal to zero. 
We will only consider totally positive sequences with cq = 1; clearly, multiplication 
of all members of a sequence by the same positive number does not affect total 
positivity. 

The following statement was independently proved by Aissen-Edrei-Schoenberg- 

Whitney in 1951 [AESW], [Ed], and by Thoma in 1964 [Th]. An excellent exposition 
of deep relations of this result to representation theory of the infinite symmetric 
group can be found in Kerov's book [Kej. 

Theorem 4.2. A sequence {c„}^o» = 1, is totally positive if and only if its 
generating series has the form 



n=0 



, ni>i(i+A») 
ni>i(i-Q!iu) 



=: TP, 



(u) 



for certain nonnegative parameters {aj}, {Pi} and 7 such that X)j(aj + /3j) < 00. 

Equivalently, an algebra homomorphism ^ : A — > C takes nonnegative values on 
all skew Schur functions if and only if the sequence {(f){hn)}n>o is totally positive, 
that is, 

n>0 

for a suitable choice of parameters (a,/?, 7). 
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The function TPa.fj^^{u) is meromorpliic in m, and it is holomorpliic and nonzero 
in a small enough neighborhood of the origin. In order to satisfy the convergence 
conditions, see Remark 2.3, we will actually have to use only the specializations for 
which TPa,f},j{u) is holomorphic and nonzero in a disc of radius greater than 1. 

Let us call a specialization a of the algebra of symmetric functions admissible if 
H{a;u) = TPa^ff^-yiu) for certain parameters (a,/3, 7) as above, and ai,Pi < 1 for 
aU i > 1. 

We say that a specialization a is a union of specializations aW,a(2),...,a('") if 

H{a; u) = H{a^^^;u)H{a'^'^^;u) ■ ■ ■ H {a<-"'^ ; u) . 

Note that unions of admissible specializations are admissible. We will use the 

notation a = l+J^j^ a^'^). 

Proposition 4.3. For any N > 1 and any admissible specializations a[l], . . . , 
a[N],b[N] of A, denote 

r r 

a{a,T]= [+J a[k], b{(T,T] = [+J b[k], a = a{0,N], b = b{0,N]. 

fe=cr+l fc=cr+l 

Assume that a = b. Then for any C > there exists a unique random point process 
on {!,..., N} X Z with determinantal correlation functions 

Pnijl ) 2;i; . . . ; T„, Xn) = det[/Cxi,rj {^i ~ 2^i)]i,j = l 

{here n> 1, Tj e {!,..., N}, a;^ e Z are arbitrary), and the correlation kernel 

' J_ r iH{aia,T];C-')Hibia,T];Oy' dC 

ICa,T{d) = { ' ' 

_J_ r H{a{T,ayX-')HibiT,a];0 dC 



(7 < r, 
a > T. 



Proof. Consider the periodic Schur process with t € (0, 1) and specializations 
a[l], b[l], . . . , a[N], b[N] as in the hypothesis. By Theorem 4.2, all weights of this 
measure on are nonnegative. Mapping to 2^^'---'^^^^ via 

(A(^), . . . , AW) ^ ({Af^ - i + . . . , {Af ) - i + 1}^, J 

yields a random point process on {!,..., iV} x Z' whose correlation functions were 
computed in §2. Shifting the space variable x 1-^ [lnC/lnf~^] + x + ^ and using 
Theorem 3.1(ii) we see that the correlation functions of our process converge, as 
t ^ 1—, to those given in the hypothesis above. Lemma 4.1 completes the proof. □ 

The main result of this section is the following statement. 
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Theorem 4.4. For any doubly infinite sequences {a[A'], of admissible spe- 

cializations of A and any c € (0,7r), there exists a determinantal^* point process on 
Z X Z with the correlation kernel 



' -i- / {H{a{a,r];C-')H{b{a,r];Oy' -J^, a < r, 



K^a,T{x-y) = < 



7 H{a{r,a];C-^)H{b{r,a];0 a > r, 



where the integrals are taken over positively oriented arches of the unit circle. 

Comments. 1. The equal time values of the kernel above are exactly those of the 
discrete sine kernel on Z; for any r € Z 

^ic{x-y) _ ^-ic{x-y) sm{c{x - y)) 



i^x-y+i 2ni{x — y) tt{x — y) 



In this sense the kernels /Co-,t {x — y) provide extensions of the discrete sine kernel. 

2. The point processes in questions are clearly invariant with respect to the shifts 
of the space coordinate. If one wants the processes to be invariant with respect to 
the time shifts as well, one has to take all specializations a [A;] to be the same, and 
all specializations b[k\ to be the same. 

3. Using Lemma 4.1 wc can slightly relax the conditions on the specializations 
a [A;], b[k] by allowing the parameters ai to be equal to 1. (Recall that they are 
required to be strictly less than 1 by the definition of admissible specializations.) 
Taking trivial specializations a[k] and choosing b[k] so that H{b[k];u) = (1 — u)~^ , 
we obtain the kernel 



K^a,T{x-y) 



which is the incomplete beta kernel of [ORl]. We will also see this kernel arising 
in the bulk limit of the cylindric partitions in the subsequent sections. 

4. The choice H{a[k];u) = H{b[k];u) = e™"''*'= " yields the extension of the 
discrete sine kernel obtained in [B02, Theorem 4.2]. 

Proof. It suffices to prove the statement for the processes on {1, . . . , TV} x Z and 
finite sequences of admissible specializations. Indeed, then one can 

just embed such processes in 2^^^ and take the limit as iV ^ c» using Lemma 4.1. 

In the finite N case we use Proposition 4.3 to construct a process on {1, ... , 2A^+ 
1} X Z with extra 2(7V+ 1) specializations a[N+l], b[N+ 1], . . . , a[2N+l],b[2N+ 1] 
defined by 

a[N + k]^b[k], b[N + k]=a[k], k = l,...,N, 
H{a[2N + 1]; u) = H{b[2N + 1]; u) = e^", M > 0. 



*i.e., with determinantal correlation functions. 
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Then 

2N+1 2N+1 

[+) a[k] = a = b= [+J b[k], 

k=l k=l 

H{a[2N + l]:C^)Hib[2N + 1]; C) = e^^*^^. 

Choosing the constant C in Proposition 4.3 to be C = e'^^''°^'^, for ( on the unit 
circle we obtain 



hm 



1, > cose, 



M^oo i^C{H{a;C-'^)H{b;C)) ^ I 0, 5iC<cosc, 

0, > cose, 



M™oo l + C-iif(a;C-i)i?(6;C) I 1, 5?C<cose. 

Hence, as M ^ oo we have the convergence of the correlation functions of our M- 
dependent processes restricted to {1, ... , N} x Z to the needed values, and Lemma 
4.1 completes the proof. □ 



5. Cylindric partitions 

Cylindric partitions were first introduced by I. Gessel and C. Krattenthaler in 
[GK] . We will initially follow their paper in our exposition. 

Let A and /x be two partitions. Assume that A D /i and denote the length 
(=number of nonzero parts) of A by L A plane partition of shape X/n is a planar 
array tt of integers of the form 

7r2,^2 + l ■ ■ ■ "'S.Mi + l ■ ■ ■ "'S.Ai ■ ■ ■ "'2,^2 
^^1,111 + 1 ••• '"■(./ia + l ••• + l ••• ■"■J.Ai 

such that the rows and columns are weekly decreasing: nij > 7r,j+i and TTij > 

A cylindric partition of shape X/ /.i/d can be viewed as a plane partition with an 
additional relation between the first and the last rows. This relation depends on 
an integral parameter d and it has the form ttij > TTi j+rf for all j. In other words, 
a cylindric partition has to remain a plane partition when the last row shifted by 
d to the right is placed on top of the first row. Here is an example of a cylindric 
partition of shape (8, 6, 3)/(3, l)/4. On the second picture the shifted last row (in 
bold) is placed on top of the first one. 

11 9 1 

7 5 2 1 
10 10 6 5 1 1 
11 9 1 

Note that d has to be greater or equal to /xi . 

The norm |7r| of the cylindric partition tt is the sum of its elements. In the 
example above |7r| = 68. 

For our purposes it is more convenient to parameterize cylindric partitions dif- 
ferently. Namely, let us read them along the lines with fixed content j — i (these 



7 5 2 1 
10 10 6 5 1 1 
11 9 1 
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lines are parallel to the diagonal which has content zero). On each line we observe 
an ordinary partition, and partitions on neighboring lines are different by adding 
or removing a horizontal strip. The necessary number of fixed content lines to be 
taken into account is equal to N := I + d, and the content can be considered as an 
element of the cyclic group Z/A^Z of order N. 

The example above leads to the following (periodic with period N = 7) sequence 
of partitions: 

... ^ (11, 2, 1) ^ (9, 1) ~< (10, 1) y (10) y (6) ~< (7, 5) y (5, 1) ~< (11, 2, 1) ^ . . . 

The or 'V' relation of the neighboring partitions depends on the boundary 
of the Young diagram /z: horizontal edges correspond to 'V and vertical edges 
correspond to . 

It is impossible to reconstruct A and jj, from such a sequence of partitions: in the 

example above if wc remove the last row from the second picture then the resulting 
partition sequence will be the same while we would be looking at a cylindric par- 
tition of type (6,6,6)/(3, 2)/4. However, d and I remain invariant — d is equal to 
the total number of "-<" and I is equal to the total number of 'V" in a period. 

Let us encode the sequence of 'V and by assigning to a cylindric partition 
a periodic sequence of I's and — I's; I's correspond to :^'s and — I's correspond to 
-('s. Thus, the example above produces (. . . , —1, 1, —1, 1, 1, —1, 1, —1, 1, . . .) (this 
sequence is periodic with period N = 7). We will call this sequence the profile of 
the corresponding cylindric partition. Clearly, the profile depends only on /U, d, and 
I, but not on A. 

We say that the profile is marked if there is a marked —1 in each period, and 
the distance between any marked —I's is a multiple of the period. In other words, 
if wc consider profiles as maps from Z/NZ to { — 1,1} then marking corresponds to 
choosing a distinguished element of Z/iVZ in the preimage of —1. 

In order to associate a marked profile to any cylindric partition we will mark the 
— I's corresponding to the relation of partitions in lines with content d and —l + l. 
Thus, in our example the marked —I's corresponds to (5,1) C (11,9,1) and the 
marked profile is 

C -1* 1-111-11 -1* 1 ) 

Marked profiles are in one-to-one correspondence with triples {fx, d, I) with d> jii 

and I > 

Let us associate to any marked profile two periodic sequences A[k], B[k] with 
period N and elements or 1 such that th(^ ciiffca-cncc A[k] — B[k] gives the element 
of the profile with distance k from the marked —1. In the example above 

A= (...,0*, 1,0,1, 1,0, 1,0*, 1,...), 

B = (...,r, 0,1, 0,0,1, o,r,o,...), 

where we marked >l[OmodA^] and B\f}modN]. As was mentioned above, we have 

N N 

Recall that a Young diagram (equivalently, a partition) v can be obtained from another Young 
diagram k by adding a horizontal strip (notation k v or v )- k) ii and only if Vi > Ki > Vi+i 
for all « > 1. 
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For any integer m let m{N) be the smallest positive integer such that 



m = m{N) mod N. 



For instance, l(iV) = 1 and -l(iV) = iV - 1. 

In what follows we will be interested in probability measures on cylindric parti- 
tions TT with fixed /i, d, and I (equivalently, fixed marked profile), whose weights are 
proportional to s''^', where s € (0, 1) is a parameter. We will call these probability 
measures uniform as the weights of the cylindric partitions with the same norm are 
equal. The next statement provides the partition function for the weights sI'^L 

Proposition 5.1. For any partition and integers d > fi\, I > ^(/u), the following 
identity holds: 



E E 

{\:Z{X)<1} {ttitt has shape A/ju/d} 



nr 



1 



n 



1 



Jp-q)(N) + (n-l)N 



pel.JV: A[p\ = l 
geMV:B[g] = l 



where N = d + I, and A[k] and B[k] are sequences of O's and I's associated to 
(/X, d, I) as described above. 

Comments. 1. It is convenient to view the function (p — q){N) as an array on the 
N X N torus with rows and columns parameterized by p and q: 



N 
1 
2 



N- 
N- 



N-1 

N 
1 



2 N ' 
1 N- 



N -2 

N-1 
N 



3 N ■ 
2 N- 



■4 
3 



1 1 

2 
3 



N N-1 
1 N 



In the formula above we choose rows p such that = 1 and columns q such that 
B[q] = 1. Note that the sets of rows and columns chosen in such a way form a 
disjoint splitting of {!,..., N} into two sets; the first set contains d elements and 
the second set contains / elements. Hence, the total number of factors in the second 
product equals dl. 

2. The right-hand side of the formula for the partition function depends on 
the profile of {^,d,l), but not on the marked profile. This is in agreement with 
rotational symmetry of the problem. The sum in the left-hand side of the formula 
may be viewed as the sum over all cylindric partitions with fixed marked profile. 

3. If d = then ji must be empty, and all rows of the corresponding cylindric 
partitions become identical. Accordingly, the second product in the right-hand side 
disappears, and wc recover the formula for the partition function of the weights 
g(# of boxes)-Ar ordinary partitions. 

4. If d, / » 1 and A = (1, 1, . . . , 1,0, . . . ,0*), B = (0, 0, . . . , 0, 1, 1 . . . , 1*) then 
jj. = d^ and the corresponding cylindric partitions of small enough norm arc in one- 
to-one correspondence with similar plane partitions. Accordingly, looking at the top 
right corner of the matrix in Comment 1 above, we see that the partition function 
for small powers of s looks like that for the plane partitions: Y\n>i{^~ ^")~" ■ In the 
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limit d,l ^ 00 one recovers the celebrated MacMahon's formula for the partition 
function of the weights 5'^°''™ on plane partitions. 

5. The case of largest rotational symmetry corresponds to ^ = (1, 0, 1, 0, . . . , 1, 0*) 
and B = (0, 1, 0, 1, . . . , 0, 1*). Then I = d, N ^ 2d, and = {d, d - 1, . . . ,2, 1) is 
the staircase partition. In this case the formula for the partition function slightly 
simplifies to give 



11 I _ g2nd 11 n _ g2m-l- 



n>l 



g2nd 11 l\ _ g2m-\+2{n-\)d\d ' 
m=l ^ ' 



Proof. It is well known that the skew Schur function s^/^ specialized at a single 
variable x is nonzero if and only if A ;^ /i, in which case s\/^j^{x) = x''^'"''^'. 

This observation immediately implies that cylindric partitions in question arc 
in one-to-one correspondence with trajectories of the periodic Schur process with 
period N determined by the specializations (see §1 for notations) 

a^[k] = ^A[k], bn[k] = ^B[k], k=l,...,N, 

or, in different terms, 

H(a[A;];u) = (l-s'=u)-^W, H{b[k];u) = {I - s-''u)-^^''\ k=l,...,N. 

Namely, wc read the cylindric partitions along the lines with fixed content, and the 
resulting ordinary partitions form a periodic Schur process. Recall that following 
the marking rule we denote the partition coming from the line with content (m 
mod TV) by X("^+'-\ In particular, the marked place of the profile corresponds to 
the relation A(°) -< A^^). 

The formula of Proposition 5.1 is then a direct corollary of Proposition 1.1. □ 

The one-to-one correspondence of the cylindric partitions with a given profile 

and trajectories of the periodic Schur process explained in the proof above also 
implies that we can use the results of §2 to obtain the correlation functions of the 
uniform measure on cylindric partitions. 

We define the correlation functions of the uniform measure on cylindric partitions 
with a given profile as the dynamical correlation functions of the corresponding 
periodic Schur process. 

Proposition 5.2. The correlation functions of the uniform measure on cylindric 
partitions with a given profile {A[k], B[k]}^^i have the following form: For any 
n > 1, Ti, . . . ,r„ e {1, . . .,N}, xi,...,Xn€ U , 



.„(n,..;...r„,.„)=/ det[i.(..,..;.„.,)lr,=J 



(z; i) dz 



where t = s^ , 
K{a, x; T, y) = 



n„>i(l-i")' / / F(a,C) ^3«r;;i) dCdr^ 



'3(z;t)(27ri)2 JJ^F{t,ji-^) e^{-Cnt-^-t) C+^r^v+h 
the integrals are taken over circles centered at the origin with radii satisfying 



s'^+i < Id < s", s-^ < Ir?! < s-^-\ |Cr?| > I if a 



r, 
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and 

/ r N 

F(r,C)=exp^— — ^K^B[fc](.-'=Cr+ E ^^^(^"^-'^0" 

n>l ^ ' \fe=l k=T+\ 

T N \ 

k=l k=T+l J 

Proof. The statement follows from Theorem 2.2, Proposition 2.1, and Remark 2.4. 
The only new effect that we have to be careful about is the choice of the contours of 
integration, because the function F{t, () is correctly defined by the formula above 
only in the ring .s^+^ < |^| < s"^, which explains our conditions on the contours. □ 

Remark 5.3. The integration over z in the formula for the correlation functions 
above can be carried out explicitly to produce a multivariate integral formula similar 
to that in Corollary 2.8. 

Remeirk 5.4. The function F{t, () can be written in terms of infinite products: 

_n i^'^VCt)^ J (.VC;t)^ 

^, kei,T:A[k] = l keT+l,N:A[k] = l 

^(t, C) = 



fcel,T:B[fc] = l keT+l,N:B[k] = l 

where, as usual, (a;t)oo = nn>i(-'- ~ at"). This formula also provides the analytic 
continuation of F{t,() from the ring 5^"+^ < |C| < s^- 

6. Bulk of large cylindric partitions 
with finite or slowly growing period 

In this section we compute the local limit of the correlation functions of the 

uniform measure on cylindric partitions with fixed profile {^[fc]. -B[fc]}^-^ in two 
cases: when N is fixed and the profile does not depend on the small parameter 
r = \nt~^ = N\ns~^, and when A'' is growing in such a way that r still tends to 
as s 1 — . 

In both cases an important role is played by the parameter 

^_ I _j:k=iB[k] 



which we call the slope of the profile. In the case of growing TV in order to have a 
limit we will need to assume that the slope tends to a limiting value strictly between 
and oo. 

We need to do some preliminary work before stating the results. 

For any >c > 0, let T^^ be the contour in the complex plane defined by 



sin(l+>r) 
sin ^cf> 



sin T±z 

- . e'ny>e(-7r,7r) 
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This is a piecewise smooth closed curve which has a corner-like singularity at 
z = 1 as (j) ^ -'-TT>t (p ^ ±7r. The origin is located inside P;^. The curve is also 
symmetric with respect to the real axis, it intersects the negative semiaxis at the 
point z = —x^^, and T^-^ is located inside Tj^^ if > ^2- Note that Fi is the unit 
circle, which is the only case when the curve is smooth at 2 = 1. Whenever we use 
as the integration contour, we will assume that it is oriented counterclockwise. 

Our interest in the family {T^} is explained by the following statement. We use 
the principal branch of the logarithm function below. 

Lemma 6.1. For any k > Q, set ^^{z) = >fln(l — z) + ln(l — z~^). This is a 
holomorphic function on C\ M>o, and 

{z e C \ R>o I fU^) e M} = K<o u . 

Further, ^fi({z) < if and only if z is located inside T ^ and below the real axis or 
outside and above the real axis. 

On the curve the function f^ equals 

fh - si£i2±iM e'^) = xln ( ESi^liM) + In f EI^LiLtlR) . 

J \ sin >c<p J \ sin >f(p / y sin (p J 

This function increases on (—3;^, 0). decreases on (0, ^jqrx); md its maximal value 
on Tit equals 

/^(-x-i) = ln(l + x) + xln(l + x-i). 

Proof. We have 

df^jre^'^) _i/ l + x\ _i/ (l + xKl-re-^X 

= r X — — = r \ K ^-7: , 

dr V V \l-re"f\^ J 

whence on any ray {z = re*"^ | r > 0}, the function '^f^{z) is strictly decreasing as 
a function of r for (f £ (0, tt), and it is strictly increasing for (p G (— tt, 0). Thus, on 
any such ray there cannot be more than one point z such that f:n{z) G M. 
On the other hand, for ^ = 1 — '^'"^^^^ e"'^ e we compute 

~ _ sintp i(K+l)4> 1 _ smxcf) -i(x+l)<t> l ^-1 _ sin(^+l)0 ^-i>«^ 

sin >C(p ' sin ' sin ' 

which implies the first statement of the lemma and the formula for f^ on the curve. 
Since 

f siri(j£+l)0y _ sin(2jf+l)(/>- {2>c+l) sin (p „ 
sin>e4, J ~ 2(sin3<(/i)2 ^ 

(sin(>f+l)0\' >f sin(2>f-t-l)(/>— (2^tf-|-l) sin , p. 

sin"^ J ~ 2(sin"^p ^ 

for (/) G (0, -^yi)-, and f^ restricted to is an even function of 0, the last claim of 
the lemma also follows. □ 

For any profile denote 

T T 
fc=<7+l fc=0-+l 

Then d = A(0, TV], / = B(0, TV], and x = B(0, iV]M(0, TV]. 
Now we can state the main result of this section. 
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Theorem 6.2. In the limit s ^ 1—, the correlation functions of the uniform 
measure on cylindric partitions with a given profile {A[k],B[k]}^^-^ have a limit 
in the following sense: Choose . . . ,a;„(s) € Z + ^ such that as s — > 1— , 

rxk{s) — > 7 for allk = 1, . . . , n and som,e 7 S M, and all pairwise distances Xi—Xj = 
Xi{s) — Xj{s) remain constant. Then for any ri, . . . , r„ e {1, . . . , N} 

\hn_ pn{Ti,xi (s); . . . , r„, a;„(s)) = det [JCp^l^. {xi - Xj)] "^.^^ , 



where the correlation kernel has the form 



dC 



2ni §T, 



l + e-''(l-C)"'(l-C"^)^'' C^-f+i 



a < T, 



a > T. 



Comments. 1. By Lemma 6.1, the expression (1 — C)'(l ~ C ^)'^ is nonnegative on 

r^, which shows that the kernel /Cafr is correctly defined for a < t. On the other 
hand, for u > t we can rewrite the formula for the kernel as 

1 / e-'il-Cy~^^^''"H^-C~y~^^-''^ dC 



l + e^(l-C)'(l-C^)'' C^-y+i 



and this integral also always makes sense. 
2. The global limit density function 



1 dQ 

+ e7(l-C)'(l-C-^)'^T 



docs not depend on r, which means that it is invariant with respect to rotations 
of the cylindric partitions. Note also that /9(7) depends on the fixed profile only 
through d — A{{),N] and I = 5(0, iV], or, in other words, through the period 
N = d + I and the slope x = l/d. 

3. The simplest nontrivial (i.e., not coinciding with the uniform measure on ordi- 
nary partitions) example is d = I = 1. Then V ^ is the unit circle, and it is not hard 
to evaluate the global limit density function and the corresponding hypothetical 
limit shape (see Comment 3 after Theorem 3.1 for explanations) explicitly: 

0(7) = ^ , v(u) = u + A arctanh(Vl + 4e") + ^ . 
y/1 + 4eT ^ 

For larger values of d and I explicit integration is also possible but as d and I grows 

it becomes increasingly tedious. 

4. A formal application of Theorem 3.1 to the uniform measure on cylindric 
partitions produces the correct integrand in the formula for the limit kernel, but 
the integration contour is different unless x = \. For x ^ \ using the unit circle 
(as Theorem 3.1 suggests) instead of the curve may lead to a wrong answer! 
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Theorem 6.3. In the limit 

s — > 1— , d^oo, 1^00, r = — A/'lns — > 0, i/rf— >x>0, 

the correlation functions of the uniform measure on cylindric partitions with a given 
profile {A[k], have a limit in the following sense: Choose xi{s), . . . , x„(s) e 

Z+ i such thai as s 1—, rxk{s)/N —> 7 for all k = 1, . . . ,n and som,e 7 G R, and 
all pairwise distances Xi—Xj = Xi{s) — Xj{s) remain constant. Also choose the time 
moments ti(s), . . . ,t„(s) such that as s — > 1 the pairwise distances Ti{s) — Tj{s) 
remain uniformly bounded. Then 

Pn{n{s),Xiisy, . . . ,Tn{s),Xnis)) = det [/C^J^^^^^.^^^ (sJi - X^)] "^.^^ + o(l), 

where the estimate is uniform over any set of profiles with the slope uniformly 
convergent to >c as s ^ 1, and the correlation kernel has the following form: 

- y) = S.-y,o ^f 7 < 7o(x) := - , 

and for 7 > 70 (>f), using the notation (^{(f>) = 1 — ^^'^{^+^^'1' g«0^ }iave 



C(0) '> 
C(0) ^ 



where both integration contours leave the origin on their left sides, and the constant 
(j) = ^(7) g (0, from the limits of the integrals above is uniquely determined 
by the relation 

-(1 + x)7 = xln ( SilliilliM) + In ( sMil^) . 

\ ' / I \^ sin>t0 J y sin0 J 

Comments. 1. The existence of the unique 0(7) for 7 > 70, as weU as nonexistence 
of (j) satisfying the relation above for 7 < 70, follows directly from Lemma 6.1. 
Using this lemma it is also easy to sec that the statement of Theorem 6.3 can be 
formally deduced from Theorem 6.2. Indeed, 

ln(e^(l - C)'(l - C'Y) =N{j^ + ^ ln(l - C) + ln(l - C')) 

= T^(^7 + /4C))- 
Thus, for Q G T^, A'' — > 00, 7 — > 00, and 'y/N — > 7 we have 

JV^oo [ +00, /^(C) > -(1 + x)7. 



Substituting this limit relation in the formulas of Theorem 6.2 we recover the for- 
mulas of Theorem 6.3. 
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2. The global limit density function for 7 < 70 (x) is identically equal to 1, while 
for 7 > 7o(x) it is equal to 



C(0) 



C 



Note that this expression is independent of t, which reflects the rotational invari- 
ance of the limit density. As 7 ^ 70 we have 0(7) and p(j) 1. On the other 
hand, as 7 — > +00 we have ^{j) — > and ^(7) — > 0. 

The fact that /9(7) = 1 for 7 < 70 means that the random Young diagrams in 
question have the lower edge — the event of having columns of length substantially 
greater than iV7o(x)/r has vanishing probability. 

3. For X = 1 the formulas simplify, and we obtain 

-In i^W^) + 1- (^^^^) = 21n(2cos0), 70 = - ln2, 

( ^ = / f arcsin(%l), 7 > - ln2, 
\ 1, 7<-ln2. 

Since |:arcsin(^^) = i arccos(l — ^-j— )> in this case the global limit density 
coincides (up to rescaling of 7 by 2) with that for the largest section of random 
plane partitions with imiform measure, sec §3.1.10 of [ORl] with r = 0. 

4. Equal time values of the limit correlation kernel equal 

CW 

^ ( _ f dC _ /sin,.0A"-^ ^in((7r-(l + x)<A)(ar-j/)) 

^r,Ax y) - J -[sin4>) ^(a- _ 

C(¥>) 

The prefactor ( ^^^1 can be ignored as it cancels out in the determinantal 
expression for the correlation functions, and we obtain the discrete sine kernel. 
Hence, the full kernel K.ct,t{x — y) can be viewed as an extension of the discrete 
sine kernel. Let us show that this extension belongs to the family of extensions 
constructed in Theorem 4.4 above. For cr < r we have 

W) 

(1 - o^-'^rfC 



/C<,^(a;-y) = (-l)^(<^'^l-^ / ^ 



V+1+A((7,t\ 



while for (7 > T 



/C.,(a.-j/) = -(-l)^(-<^l(^) — j 



Qx-y+l+A{a 



t((l + >«)<<.-7r) 



(^x-v+1-A(t,<t] 
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Once again, the prefactors can be ignored, and after shifting the space variable x 
at time r by a; i— > x + A{T, r] for a fixed T and all r, we see that the correlation 
functions are determinants of the kernel 



i(,r-(l + «),;,) 



V sm >ccb ^ 



2Tvi J i-x-y+\ ' O ^T, 

ei((i+«)*-^r) S 



2¥i J Ax-y+l ' cr > T, 



where the integrals are taken over positively oriented arches of the unit circle. 

For X = 1 this is exactly the incomplete beta kernel arising in the bulk limit of 
uniform measures on plane partitions, sec [ORl]. 

If X > 1 then -^-^ < 1 and the kernel coincide with one of the stationary 
kernels afforded by Theorem 4.4; one has to take H{a[k];u) = 1, H{b[k];u) = 
(1- 4H^ru)-i for all fceZ. 

V sin >f(p ' 

Finally, if >f < 1 then ^^'^^^ > 1, and one more shift of the space variable x 
at time rhy xi-^ x — T + t (or, equivalently, the shift a; i— » x — B{T, r] of the 
initial space variable x at time r) for a fixed T and all r, brings our kernel, up to 
irrelevant prefactors, to the kernel of Theorem 4.4 with H{a[k];u) = (1 — u)~^ 
and H{b[k];u) = 1 for all k€Z. 

Observe that all these kernels are invariant with respect to time shifts. Thus, 
the bulk behavior of our cylindric partitions with slowly growing period N is rota- 
tionally invariant and independent of the initial profile except for its slope, after we 
perform the above shifting of the space variables. However, these shifts do depend 
on the profile. 

Proof of Theorems 6.2 and 6.3. As in the proof of Theorem 3.1, let us introduce a 
new integration variable ^ = Cv and rewrite the kernel from Proposition 5.2 in the 
form 

. ^_ n„>i(l / F{aX) 03{z^;t) dCd^ 
i\ [a, X, T, y) - 



e^{z;t){2niY J J F{T,a-^) 03(-^i-^;i) (-y+^^y+i 



where, by Remark 5.4, 

F{a, ^ Uke—:B[k]=ii^ - s^'^C/O nfcg^+uv:A[fc]=i(l " ^VO 

11 (<l-k^-r■f^.. 11 



kel.N:B[k] = l ^ ' k(i\,N:A\k\ = \ ' ' 

Here, according to Proposition 5.2, the integrals are taken over circles centered at 
the origin with radii satisfying 5^^+^ < \C\ < s'^ and s'^^^ < |CI/ICI < s'^ plus the 
additional condition |^| > 1 if cr = r. 

The proof proceeds along the same lines as that of Theorem 3.1. The only major 
obstacle is the unboundedness of the ratio F(cr, (t, CC~^) on the integration 
contours. Recall that in Theorem 3.1 we imposed the condition Ak = Bk which 
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ensured the needed boundedness when \C,\ and |^| were close to 1. Here we only 
have such "self-adjointness" in the case x = 1; for general x > we need to deform 
the ^-integration contour. 

Clearly, the first factor in the formula for F{aX)/F{TXi~^) above remains 
bounded as s — > 1 as long as the factors in the denominator do not approach zero. 
Furthermore, as ^ ^ 1 and s ^ 1, it tends to (1 - C)^''''^l(l - C"^)^^''''^^ for a <t, 
and to (1 — C) (1 ~ for a > t. As we will see, only integration 

over infinitesimally close to 1 ^'s yields a nonzero contribution, and this asymptotics 
gives the corresponding factors in the final formula for the correlation functions. 

The second factor needs more attention. Recall that the dilogarithm function is 
defined in the unit disc by the power series 



which shows that 2: = 1 is the branching point of this function. We will consider 
dilog(z) as the holomorphic function on C \ (1, +00) defined by the integral above 
with the principal branch of ln(l — a;). Note that the jump of dilog(2;) across (1, +00) 
is purely imaginary, which means that 5idilog(2:) is a continuous function on C. 
Our interest in dilog( • ) is explained by the fact that 



is, up to the factor (1 — a Riemannian sum for the integral representation 

of dilog(a;). Thus, using the mean value theorem to estimate the remainder, for 
a; e C \ [1, +00) we obtain 

ln(a;; t)oo = -(1 - i)"^ dilog(a;) + 0(dist(l, {ra; | < r < 1})"^) , t ^ 1, 

and the estimate is uniform on compact sets. (The remainder may become large as 
X approaches the cut [1, +00).) 

This asymptotic relation shows that the second factor in the formula for the 
ratio F{aX)/ F{t, C^^^)j as f ^ 1, is approximated by 



Z(dilog(C/0 - dilog(C)) + rf(dilog(l/C) - dilog(^/C)) ~ In ^{l ln(l - C) + d ln(l - C"^)), 



which coincides, up to multiplication by dln^, with function f^tiC) introduced in 
Lemma 6.1 above. 




Its analytic continuation is provided by the integral representation 






Note also that, since dilog'(a;) = — ln(l — x)/x, as ^ ^ 1 we have 
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Lemma 6.4. For any x > and ^ 7^ 1 on the unit circle, set 

/.,aC) = ^(dilog(C/0 - dilog(C)) + dilog(l/C) - dilog(C/C). 

Then f^,^ is a holomorphic function on C\{M+U^M+}, its real part is continuous 

on <C, and 

{C e c 1 3?/^,j(C) = 0} = • R u r^,j, 

where ,c is a piecewise smooth closed curve which encloses the origin and inter- 
sects each ray e*'''M+ ^ • M a single point. Further, ^f>c,^{C) < if and only 
if ( is inside F^,^ and to the left of the line -R, or Q is outside T^ ^^ and to the 
right of the line • M. 

The curve Fj^^^^ lies inside F^^^^j if > X2: the curve Fi^^ is the unit circle. 
Also, as ^ 1, we have T^^^ in the sense that the intersection points ofT^^^^ 

with rays e"^]R+ converge to the corresponding intersection points ofT^. 

Proof. The fact that takes real values on • ^ follows from the relation 
dilog(a;) = dilog(S). On the other hand, on any ray we compute 



dr 



This expression is negative if ( is to the right of the line • R, and it is positive 
if C is to the left of the line y/i ■ R. 
Using the asymptotic relation 

lim Sfdilogf-re'^^M - dilog(-re*'^=)) = - ^l), -n < V'l, V'2 < tt, 

1 >+(X) 

we SCO that the limits of fx,^{re^'^) as r ^ 0+ and r +00 exist, and if e''^]R+ ^ 
^/^ ■ R then they are nonzero and have different signs. Since the real part of 
the derivative along any such ray is sign definite, this implies the existence and 
uniqueness of the needed intersection points as well as the inequalities on ^f^^^{Q. 

Similar arguments show that vanishes at exactly two points located in 

• ]R+ and -^/f • M_. These are the points where Fj^^j intersects the critical line 
v^ M. Thus, the curve is closed, and a standard implicit function theorem argument 
implies the smoothness. 

Since the limit values limj.^o+ fK.^{re^^) do not depend on x, while the absolute 

value of the derivative |3? '^^"'^^^.^ — ^| is an increasing (and linear) function of n, we 
see that the larger x the sooner f,t,^{re^'^) reaches as r increases. This implies 
the inclusion property of the curves r;^;,^. Finally, the convergence of r;^;,^ to T ^ 
follows from the asymptotic relation /ji,j(C) ~ ln^/>((C) as ^ ^ 1, Lemma 6.1, and 
the fact that In^ is purely imaginary as |^| = 1. □ 

Prom this moment the proof very much resembles that of Theorem 3.1. For 
^ bounded away from the point 1 we use Lemma 6.4 to deform the ^-integration 
contour to the one where the ratio F{a,C,) / F{t,QS^~^) remains bounded as i ^ 1. 
The result of Proposition 3.2(i) then implies that the integral over such ^ tends to 
zero. On the other hand, for ^ close to 1 we deform the ^-contour to T^t and use 
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Proposition 3.2(i) together with the asymptotics of F{a,Q/F{T,(^ ^) to see that 
(cf. the formula before Lemma 3.3) 



K{a,x;T,y) 



Hkel,r:B[k] = l(^ Ilkea+l,N:A[k] = li^ ^/C) 



'^(i^r Jt^ n.ei;^:B[.]=i(i - c) nfee7+T;iv:Aw=i(i - VC) 

r+ooife g^(lnz-7-nn(l-C)-ciln(l-C-l)) ^^^^ 
f,u _ p-u (-x-y+1 



r 

J — ( 



as t ^ 1, e > is small enough, the sign in ztie is "+" for a < t and "— " for cr > r, 
and 7 is replaced by in the case of Theorem 6.3. The integration variable u is 
related to (, via u = —i'jTr~^ In^. 

The integral over u is evaluated by Lemma 3.3 to give 

±in 



l + {z-^e^{l-Cy{l-C-')df^ 



Finally, similarly to the proof of Theorem 3.1, the integration over z (see the formula 
of Proposition 5.2) is shown to be asymptotically equivalent to substituting z = 1. 

In the case of Theorem 6.3 the argument is exactly the same except for the fact 
that 6^(1 — C)Hl ~ C""'^)'^ converges to either or +oo, sec Comment 1 to Theorem 
6.3. However, Lemma 3.3 can still be applied because its estimates remain uniform 
as long as (in the notation of Lemma 3.3) e5fa is uniformly bounded and |9a| is 
bounded away from 1. Since e can be made arbitrarily small, 3?a is allowed to 
converge to ±oo. 

We skipped the discussion of several technical issues here. One should make 
sure that the deformation of C-contours happens in such a way that the poles of 
the ratio F{a,Q/F{T,C£.^^) Stic not passed in the process; the remainder in the 
approximation of ln(a;, t)oo by dilog(x) may grow as x approaches cut x € (1, +oo), 
and one needs to control the growth of the derivative of the remainder term in order 
to apply Lemma 3.3. All these issues can be resolved. However, the arguments are 
rather tedious although fairly straightforward, and we omit them. □ 

7. Bulk of large cylindric partitions 
with period of intermediate growth 

In the previous section we have seen that if the period A'' is small comparing 
to |lns|~'^ then the global density function for the uniform measure on cylindric 
partitions with a fixed profile, when the partitions are scaled by | lns|, converges to 
a rotationally invariant limit which depends only on the average slope of the profile. 

In this section wc consider the case when t — has a nontrivial limit strictly 
between and 1. Equivalently, N In ^ \nt~^ G (0, +oo). This situation is much 
more complex because, as we will see, the global density is no longer rotationally 
invariant. 

It is worth noting that if N grows fast enough (so that t ^ 0) then the random 
cylindric partitions split into disjoint independent random plane partitions located 
at the "corners" of the profile provided that these corners are deep enough. 

We will consider in detail two cases: the most rotationally symmetric one with 
and being periodic of finite period (in the notation of §5 this means 

that fj, is the staircase- like partition), and the most rotationally asymmetric one 
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with {j4[fc]} and _B{[fc]} consisting of two blocks of zeroes and ones (this means 
that n is either empty or it is the rectangle dx I). 

In both cases we will assume that the average slope >c is equal to 1. This will 
allow us to obtain asymptotic results without substantially deforming the contours 
in our integral representation for the correlation kernel. We hope to consider more 
general cases in a subsequent publication. 

Proposition 7.1. Consider the uniform measure on cylindric partitions with a 

given profile {A[k], B[k]}^^-^ and assume that the sequences {A[k]} and {B[k]} are 
periodic with an even period M and A{0, M] = B{0, M] = M/2. Then as 

s^l-, N^oo, s^^ie(0,l), 

the limit of the correlation functions of this measure is described by Theorem 6.3 
above with x = 1 (see also Comment 3 after Theorem 6.3). 

Proof. For the sake of convenience let us consider the case when the large period 
AT is a multiple of the small period M; N/M e Z. This assumption is by no means 
necessary and can be easily removed. 

Also, using the rotational invariance, we may assume without loss of general- 
ity that all the time moments ti(s), . . . , r„(.s) from the statement of Theorem 6.3 
remain positive and uniformly bounded as s ^ 1. 

As in the proof of Theorems 6.2 and 6.3, we start with the integral representation 
for the correlation kernel K{a,x; r, y) given in Proposition 5.2. We have 

F{a,C) ^ nfceT;7:B[fc]=i(l " ^'V^?) IlkeT:^:A[k]=ii^ ' ^S) 

nr,v-') n.ei;^:sw=i(i - «-'c) nfeeT;^:A[fe]=i(i - ^vo 

kel,N:B[k]=l kel,N:A[k]=l 

Using the periodicity assumption, we can rewrite the second factor as 
^1 (gM-fc^ M) _L1 {skrj-s^U 

fcel,M:B[fe] = l feel,M:A[fe] = l 

(here the base t of all products {x;t)oo changed from t to s^). 

Starting from this moment the proof almost literally repeats the proof of Theo- 
rem 2 in [ORl], see §§3.1.3-3.1.9. Approximating {x; s^') by the dilogarithm func- 
tion (see the previous section) and keeping in mind that A{0, M] = B{0, M] = M/2, 
we obtain that as ,s ^ 1 the second factor in the formula for the ratio F{a, Q /F{t, rf) 
above is approximated by 



exp 



M(dilog(C) - dilog(l/C) + dilog(?y) - dilog(l/r/)) 
2(1 - s^) 



Note that 1 — s^ can be replaced by M\ lns| as s — > 1. 

The only other two factors in the integral representation of K{(7,x;T,y) which 
may not have a finite limit as s ^ 1 are C,~^~^r]~y~^ . Since x,y ~ 7/|lns|, we 
can include these factors into the exponential above to obtain 



exp 



'dilog(C) - dilog(l/C) - 27lnC + dilog(77) - dilog(l/77) - 27 In 77' 
. 2|lns| , 
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Following §3.1 of [ORl] we see that the function 

5i(dilog(C)-dilog(l/C)-27lnC), 

which clearly vanishes on the unit circle |C| = 1, is negative slightly inside the unit 
circle if — ^ arcsin(^^) < arg^ < ^ arcsin(^^) and shghtly outside the unit circle 

if |argC| > ^arcsin(s^); the points arg( ~ ±i arcsin(S-2~) ^ith |C| = 1 are the 
critical points of the function dilog(l/C) — dilog(C) — 27lnC. 

We would like to deform both C,- and ry-integration contours to put them in the 
domain where the exponential above converges to zero, i.e. where 

3?(dilog(C) - dilog(l/C) - 27lnC) < 0, 3?(dilog(7?) - dilog(l/7?) - 27 In 77) < 

(this corresponds to using the contour 7< from [ORl]). However, there is an ob- 
stacle: the hmction Q'^{-~C,r\t^^: i}j, which enters the denominator of the integrand 
in the formula for K{a,x;T,y), vanishes when (^rj = 1. Thus, while deforming the 
contours we have to add the residues corresponding to the poles (rj = 1. 

Consider the case cr < r first. Then Proposition 5.2 states that the C,- and 
?7-contours are such that ry-contour contains the l/(^-contour. Since our desired 
contours go inside the unit circle when the arguments of C and r] are inside the 
interval (— ^ arcsin(^^), i a,TCsm{^^j^)) , this is where the residues have to be added. 
For 7 < — In 2 the residues have to be taken along the whole integration contour. 

Observe that 

n„>i(l-i")'^3«?7;i)F(a,C) 1 

- Res^=i/f 



Y\.kea+l,T:B[k] = li^ ~ ^ ''0 YikeT+l,r:A[k]=l ~ ^'^ / 



which asymptotically equals (1 - ()B{<y,r]^i ^-i-^A(a,T] ^^x-y+i as s ^ 1. (Recall 
that we assumed in the beginning of the proof that the time moments cr, r remain 
bounded as s ^ 1.) Thus, the residue contribution is asymptotically equal to 



27ri X 



^■^ arcsin(^-2 — ) 



^ arcsin(£-y^) ^ 



x-y+1 



On the other hand, the integral when ( and rj range over the desired contour 
converges to zero as s ^ 1, see §3.1.1 of [ORl] for additional explanations. 
The arguments for cr > r are very similar. 

Note that all the estimates are uniform in the parameter z, \z\ = 1, and the final 
asymptotic expressions are independent of z. This means that the integration over 
z in the formula of Proposition 5.2 in the limit s ^ 1 can be simply removed. □ 

Now let us consider the non-symmetric case. Take d = I = N/2 and 
{^[fc]}feLi = (1, 1, . . . , 1, 0, . . . , 0, 0), {S[A;]}f=i = (0, 0, . . . , 0, 1, . . . , 1, 1) 

V ' V ' ^ V ' V ' 

N/2 N/2 N/2 N/2 
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(we are assuming that N is even). This corresponds to the square partition fj, = 
(■^)^, and the random sequences of ordinary partitions of the form 

We will describe the limit of the correlation functions near the largest partition 
Ad) and near the smallest partition A*^"+^). (Recall that the time variable of the 
periodic Schur process or, equivalently, the index of the Young diagrams above, 
may be viewed as an element of Z/7VZ.) 

Theorem 7.2. In the limit 

s^l-, N^oo, s^^te(0, 1) 

the correlation functions of the uniform measure on cylindric partitions with the 
profile as above have limits in the following sense: Choose xi{s), . . . ,Xn{s) G Z + ^ 
such that as s ^ 1—, | lns| • Xk{s) 7 for all k = 1, . . . , n and some 7 e R, and 
all pairwise distances Xi — Xj = Xi{s) — Xj{s) remain constant. 

(i) Choose the time moments ti(s), . . . , t„(s) € {— "j, • • • , -j} such that as s ^ 1 
the absolute values \Ti{s)\ remain uniformly bounded. Then 

Pn{Tl{s),Xl{s)] . . .,Tn{s),Xn{s)) = det (oJi - 2;^)]^^^ + o(l), 

where the correlation kernel has the following form: 

4^ (a; -y) = 5{x -y) if ^ < 7r"(t) := 2 In , 
anrf/or7 >7r"(i) 



g-to s 



^uhere the constant c = 0(7) e (0,7r) from the limits of integration is uniquely 
determined by the relation 



7 = 21n 



(e-;t)c 



(a) Choose the time moments ti(s), . . . , r„(,s) e {1, . . . , iV} such that as s ^ 1 
the absolute values \Ti{s) — remain uniformly bounded. Then the asymptotics 
of the correlation functions has the determinantal form as in (i) above with the 
correlation kernel given by 



K^%x-y)=< 



S{x-y), if 7<72™(t):=21n ^ 



0, if -i>^r^{t):=2\n 



{-Vi;t)^ 



(\/i;Oc 
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and for 7 € i'^™™{t),^^^^{t)) it is given by the formula of (i) above, where the 
constant c = 0(7) € (0, tt) from the limits of integration is uniquely determined by 
the relation 



7 = 21n 



{e^^Vt-^t) 

Comments. 1. For any t G (0, 1) the functions 

{e''Vi;t), 



7i(t,c) =21n 



(e'^;t)c 



72(t,c) =21n 



(e'^t;t)c 



(e'^\/t;t)c 



as functions of c are strictly decreasing on (0, tt), and 

lim7i(t,c) = +00, lim 7i(<,c) = jf^it), 

C^O C^TT 

lim 72 {t, c) = 'yr^'it) , lim 72 (t, c) = jf^'it) . 

c— ^0 c— ^7r 

Indeed, this follows from the fact that for any a,f3 G (0, 1) and c € (0, tt), we have 



d_ 

dc 



1 - af3e' 



1 — ae* 



2a(l - g^/?) sin(c) 

11 - ae*=|4 



<0. 



The decay of 71 (f, c) and 72 {t, c) guarantees the existence of the unique c € (0, tt) 
satisfying the needed relation. 

2. The global density function in both (i) and (ii) is given by 



p(7) = < 



r 1, 7 < 75^^(0, 

^c(7), 7iT(*) ^ 7 < 7f:i"W, 



V 0, 



7 > 7-|-(t), 



where 7f"''^(t) = +qo. In particular, in (i) we see that the random Young diagrams 
have only the lower edge, while in (ii) there are both the lower and the upper edges. 

3. The functions 7i,2(i,c) can be written in terms of the Jacobi elliptic sine 
function: 



(e-;i)c 



ts 



1 



2Vk sn{^,k) sm{{ 



{e'H-t), 



(e-V<;t)c 



K=iel{Q;t). 



2ti 



sm 



Note that as f — > 1, the elliptic sine tends to 1, and the formulas of both (i) and 

(ii) degenerate to those of Comment 3 after Theorem 6.3. This agrees with the fact 
that if N is growing slowly then the global limit density function is rotationally 
invariant. 

On the other hand, as t ^ in the setting of (i) we have 



\{e'<^^Vt;tU\^ 
|(e-W;i)ooP 



1 



1 



\l-ei^h)\^^ 2(l-cosc(7)) 



or 0(7) — * arccos(l — ^2~) 
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which is exactly the behavior of the largest section of the random plane partitions, 
cf. §3.1.10 of [ORl] with r = 0. This agrees with the fact that if N grows too fast 
then the random cylindric partitions split into disjoint independent plane partitions 
located in deep corners of the profile (in this case we have only one such corner). 

Proof. The proof is very similar to that of Proposition 7.1 and [ORl, Theorem 2]. 
The only difference from the arguments used for Proposition 7.1 is in the specifics 
of the functions F{t, Q. Using Remark 5.4, for bounded r > we have 

N/2 

r ,n(.svc;*)^ 
n (s'ct)^ 

fc=0 

The first factor has the finite limit (1 — C~^)~'^) while for the second factor we have 

^/2-l _^ / ,2. 



n (^'-;*)c 



fe=0 



{xVi; s)ooixtVi; s)oo{xt'^Vi; s) 
1 " 



exp — ^ (dilog(a;i™\/t) - dilog(a;t'")) , 



m=0 



Thus, in F(t, Q we have the factor 
1 



TT-n E (dilog(Ci") - dilog(Ci™Vi) - dilog(t-/C) + dilog(t'" Vi/C)) 
in s I ^-^ 



m=0 



times the part which has a finite limit as s — » 1. 

The real part of the above sum is identically equal to zero on the unit circle 
I CI = 1. Furthermore, the derivative of this sum with respect to C, multiplied by C, 
equals (for \(\ =1) 



-21nn 



m>0 



1 - Ci" 



1 - Cf^Vi 



= 21n 



(C;t)c 



From Comment 1 above we know that for any 7 e (7™'", +00) this expression is 
strictly greater than 7 on the arch arg^ G (—0(7), 0(7)), where c = 0(7) satisfies 
the relation in the hypothesis (i) above, and it is strictly less than 7 on the com- 
plementary arch. This means that if we deform the C-integration contour so that it 
goes slightly inside the unit circle for argC e (—0(7), 0(7)), and slightly outside the 
unit circle for jarg^j > 0(7), then the factor F(t, C)C~^~^ with x 7/|lns| will 
exponentially decay on this contour. 

Similarly, for r = AT — r with nonnegative and bounded f we have 

N/2 

r nits' let) ^ 

n [sHM^ 

fe=0 
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which has the same asymptotically nontrivial part as before up to the change C — » 
t(. This means that in this case the desired ^-integration contour is obtained from 
the one in the previous case by multiplication by t. 

The proof of (i) proceeds in the same way as the proof of Proposition 7. 1 : The 
deformation of contours to the desired ones meets an obstacle of the theta function 
93{—(rit~2]t) in the denominator of the integrand vanishing for (^rj = 1 or (ij = t; 
evaluating the corresponding residue yields the limit correlation kernel, while the 
remaining integrals over the constructed contours tend to zero as s ^ 1. The 
estimates are uniform in z, \z\ = 1, and the limiting kernel is independent of z; thus, 
the integration over z in the formula for the correlation functions in Proposition 
5.2 can be removed. 

In order to prove (ii) we need to consider F{t,() where r = N/2 + t with 
bounded r. The asymptotically nontrivial part is independent of the sign of r, so 
let us take r > to be concrete. We have 

N/2 

n(i.sVC;t)^ 

fc=0 

The first factor has a finite limit as s 1, while the second factor (the ratio) 
produces 

exp -r^ £ (diiogia'") - diiog(a"' Vi) - diiog(r+vc) + diiog(i™+ix/i/c)). 



m=0 



The real part of this sum vanishes on the circle \(\ = \ft. Further, on this circle 
the derivative of the sum multiplied by C, equals 



2 In 



(Coi;i)c 



(Co\/t;i)c 



c _ C 



The remaining part of the proof is just as in the case (i) considered above. It is 
worth noting that for 7 > 7™^^ the deformation of the contours to the domain 
of exponential decay of the integrand docs not meet any obstacles, and thus the 
limiting kernel is identically equal to zero. □ 



8. On a measure of Nekrasov and Okounkov 

A large part of the material presented in this section is the result of joint dis- 
cussions of the author and Grigori Olshanski. 

Our goal in this final section is to demonstrate how the techniques developed in 
the first three sections apply to a remarkable measure on partitions introduced by 
Nekrasov and Okounkov in [NO] . 

Fix /i g C and f G (0,1). The object of interest is a (generally speaking, complex) 
measure on the set of all partitions given by the formula, see [NO, §6.2], 

M,,(A) = n (1 - n-'^ • i'^' n ^^^?70T^ ' ^ e ^- 
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Here the product is taken over all boxes of the Young diagram A, and h{D) denotes 
the length of the hook rooted at the given box. The sum of the weights M^^t(A) over 
all partitions A is identically equal to 1, see [NO, (6.12)]. Note that this measure 
becomes a probability measure (meaning that all weights are nonnegative) if /i G zM. 

One interesting feature of measures M^_i is that they interpolate between the 
uniform measure on partitions arising at /j. = 0, and the poissonized Plancherel 
measure which is obtained by the limit transition 

Here dim A is the dimension of the irreducible representation of the symmetric 
group S'l;^!- We refer the reader to [BOO] for details and further references on the 
Plancherel measures. 

Let us denote by the specialization of the algebra of symmetric functions A 
such that 

hM = ^, n>0, or H{p^,u) = {l-u)-^. 

Here {a)k = a{a + 1) = . . . (a + A: — 1) denotes the Pochhammer symbol. 
The applicability of the periodic Schur process results follows from 

Lemma 8.1. For any k,X gY choose r > max{^(K),^(A)} and set 
ki = Ki - i + r, k = Xi - i + r, i= l,...,r. 

Then 

In particular, for k = X the last factor turns into 1 and we obtain 

h{Uf - ^2 



Comments. 1. Only finitely many terms in the sums over v above are nonzero. 

2. The formula for k = A can be easily extracted from the Fock space represen- 
tation of the measure M^^t given in [NO] . 

Proof. Since the terms of the sum vanish unless u c k and z^" C A, we can restrict 

the sum to u with £{v) < r. Denote Ui = Vi — i, i ~ 1, . . . ,r. 

The Jacobi-Trudi formula for the skew Schur functions, see [Macd, §1 (5.4)], gives 



Sk/u = det[/ifei-njrj=i, sx/u = <iet[hi,-nj]lj=i- 



PERIODIC SCHUR PROCESS AND CYLINDRIC PARTITIONS 



59 



Applying the Cauchy-Binet summation formula, we obtain 

X] s^/^{p^,)sx/AP-t') = X] det[hki-nj {Pti)]lj=i det[hu-ni {p-^)]lj=x 

i/eY ni>--->nr>0 

min{/ci } 



X hki-m{P^)hlj-m{P-ix) 
rn=0 



= det 
Using the identity 



det 



E 

m=0 



— m 

{ki-m)\ {lj-m)\ 



i{fe,;} 



^ + l)fc(-M + 1); 
(fc-m)! (Z-m)! k\l\{n + k-l) 



E 



m=0 



which can be proved by induction on minjfc, 1} with k — I fixed, and the formula 
for the Cauchy determinant 



det 



jJ' + ki — Ij 



2 — -'— 

Ul^=li^^ + - Ij) 



we obtain 



I Ij l 



i=i 

Finally, we use the fact that the set of hook lengths {/i(n)}neK can be obtained 
as the union of sets {!,..., fcj} for i = 1, . . . ,r minus the set of numbers {ki — 
kj}i<i<j<r, and similarly for A. This brings us to the desired formula. □ 

Let us now associate to the measure M^f^ ^ a measure on point configurations 
(subsets) in Z' = Z + i via the map 

The correlation functions of this measure on point configurations are given by 
Pn{xi,...,Xn)= X ^ti,tW- 

AGY:{Ai-i+i}i>iD{a;i,...,a;„} 

We will also need the shift-mixed version of M^^f For any z G C\{—t^^, —t^i , • • • } 
define the measure Mfj^^t,z on Y x Z as the product-measure 



cf. §2. We associate to this measure a measure on point configurations in Z' via 
the map 

YxZ^2^', \^ {S + \i-i+\).^^. 
Its correlation functions are given by 

pf^*(xi,...,a;„)= X M^,t,z(A,S). 

(A,S)eYxZ:{S+Ai-i+i}j>iD{xi,...,x„} 

The two sets of correlation functions are related as described in Proposition 2.1. 
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Proposition 8.2. The correlation functions of the shift-mixed measure M^^t,z have 
determinantal form: For any n = 1,2, . . . and xi, . . . , a;„ € Z', 

pt'^\xi, ...,Xn) = det[K^,tAxi,Xj)]lj^i, 

where the correlation kernel has the form 

n„>i(i-n' 



e3{z;t) (27ri)2 



. ^ TT ^ 

Here both integration contours are simple positively oriented loops going around the 
origin such that 

t<\C\<l<\v\<t-\ l<\Crj\<t-\ 

and we use the principal branch of the logarithm to define (•)'*. 

Proof. Consider the periodic Schur process with N = 1 and speciahzations a = a[l] 
and b = b[l] of A defined by 

H{a;u) = {1- au)-^, H{b;u) = {1 - au^, < a < 1. 

Then by Proposition 1.1 

AeY veY n>l 

which converges to nn>i(-'- " t")'^^~^ as a ^ 1. Since /(a) fiPn) and f{b) 
fiPfj.) for any / e A as a — > 1, we see that the weights of the periodic Schur process 
above as well as those of its shift-mixed version converge to the corresponding 
weights of Mf^^t and M^^t,z- Then the fact that the sum of the weights in all cases 
is identically equal to one, implies the convergence of the correlation functions of 
the (shift-mixed) Schur process to those of M^^t and Mf^,t.z- 

On the other hand, the analytic version of Theorem 2.2 (see also Remarks 2.3 
and 2.4) gives the formula for the correlation functions of the shift-mixed Schur 
process introduced above. Taking the limit q ^ 1 in the integral representation of 
the kernel given in Remark 2.4 completes the proof. □ 

Let us now restrict our attention to the case when the measures M^i^t and Mfi^t,z 
become probability measures (in other words, all weights are nonnegative). This 
happens when /j. = i^o with /xq € 1^, and z G K+. 

Theorem 8.3. Assume that ji = i/io with /io G K, and z > 0. Then as t ^ 1 
the correlation functions of the shift-mixed measure M^^t,z have the following limit: 
Choose xi{t), . . . , Xn{t) e Z' such that ast ^ 1, \ lnt\-Xk{t) — > 7 for allk = 1, . . . ,n 
and some 7 G R, and all pairwise distances i{t) — Xj{t) are independent 

oft. Then 

limpf ft(a;i(t), . . . ,x„(t)) = det[lC^^''^''^\xi - 
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where the correlation kernel has the following form 

^ ' 27TI J l + z-ie7(l-C)-^(l-C-i 



deZ. 



IC=i| 



Under the same assumptions the correlation functions of the measure M^^t converge 
to the limiting expression above evaluated at z = 1. 



Comments. 1. For \(\ = 1 we have 



gMo(argC-7r)^ < arg C < TT, 

gMo(argC+^)^ -7r<argC<0. 



The change of sign of /iq is equivalent to the change of variable C ^ in 
integral above, which in its turn is equivalent to transposing the correlation kernel. 
Clearly, this operation does not change the correlation functions. 

2. The global limit density function for the shift-mixed case is equal to 



1 

2^ 



1 



e^ + ze" 



27r/io eT + ze 



In 



and one has to substitute z = 1 to get the formula corresponding to the non-mixed 



measure M, 



3. In the limit fiQ —> +00, 7 —> cxi so that j/jio 7, the correlation kernel 
becomes equivalent to the discrete sine kernel, cf. Example 3.4: 



lim lC'^'''''-^'\x-y) 

fiQ — >-+oo,7 — >-oo, 



0, 



5{x - y), 



sin(^ {x - v)) 
7t{x - y) 



7 > TT, 

— TT < 7 < TT, 

7 < — TT. 



The proof of Theorem 8.3 is completely analogous to that of Theorem 3.1 and 
we omit it. 
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